“Calhoun 


Institutional Archive of the Naval Postgraduate School 





Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations 1. Thesis and Dissertation Collection, all items 


1991-12 


Orbital maintenance of endoatmospheric low 
earth-orbiting satellites 


Pauls, David D. 


Monterey, California. Naval Postgraduate School 
http://ndl.handle.net/10945/26679 


This publication is a work of the U.S. Government as defined in Title 17, United 
States Code, Section 101. Copyright protection is not available for this work in the 
United States. 


Downloaded from NPS Archive: Calhoun 


Calhoun is the Naval Postgraduate School's public access digital repository for 
| (8 D U DLEY research materials and institutional publications created by the NPS community. 
«ist sia Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 


NY KNOX appointed — and published -- scholarly author. 

ia) LIBRARY Dudley Knox Library / Naval Postgraduate School 

411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 





http://www.nps.edu/library 





















~ 
: 


tn 5, 


‘+ Jeo ? + + 
Seer .?pr 
se] Cutie WPT oi 

Rar ar) UGE 
a ? "oye t 


Fol teak het 9) opt 
on tLe aed " 
at Ad se ‘ erat) 
easy “Beans 

ve 






















ND frags! Siedet atts One" 
et: edb ABE Gey (A ins 


oes 
ae ee AM i 
da Le 
we saat Hest 8 

’ Dake tek, 7 ¢4 sheeak aMn 
fas oy Pht poets age Be Foley 
Shai hee ae a oe Siete 


We LOR ps aba oe Re iee e tae8 
a PACS 


















oe 


wee 
sos 





4 


' 
a oe om # eS we pet 
7 ee att ‘ ‘ 





ue uh U 
: ca a fe ce sun 
*h ne W Nc lave 


ot Peodae 





= iz, 


es 







































bi 
Va ore | 
Uaaiens ai," 


B30 Vlas 















car hit 


Bid nd it Pareto nr8, 
he a rere tree pried) 
Weds b.6 ag. 








vl 1 buyer ‘at 
AU Nee tre 
fe 4 hat 

‘ ie 














a ceva 95 
te aga Obs 
ea 


R srderg 
aye afr | a My 
‘ Ey", 

Rea ire! 



























4% “ 
ately Poke 
ae ahet aye¥ S syiages k hore 
bis pely ‘ sry e eae? 
To AE ge Tab aay ee, ee eats 
asi be. & wal omy hy 
fol 4 


Seale 







ira 


maes 






Rapes 


Mtoe, ry} 
r 


. 1 + b oat 
eget . 
. eel Se Ce | 4 


























Pre Tas 










he Uh is 4d ae iy 
Ag 














t.2 + ¢ 

t VY de temea 

, Sead Wi NS ase ks 
a yt eet Car br ot My lad X 
4h a 5 g Fhe HOLE e 











ashy 
erat Gh 
Re fra ane thus He & LB ee i! cits ey ares Dae ee 

A BS Spe es aia’ ays A ult 4a Se acu Gs, iit a3. ae ih, : : ry 

o al ae + at”; Mates Deiat apt tere PY bs 2 F 
Ll hd aeviesthes * fees i fi. & ty or sins, eas ‘ Fedth gaghene HN eee Mages a 


2 + ste” AMER © hae dp Be aurah 
he tbe Sat eats, V5.8 le on Vga _ i ae 4 ata & fn 
art a ay a aftt, Ree Poe | 4 i pas rarer rer mimcetsl ac kh rare nen sis a eas ae ee: 


ii Neher, ste We Fake ky fad a cpg ate. fa MA Obe oe 
tay a 3 . + ay U Ti a * 
= fit aet dl roe 4 ve, eit ce atvas Py are FASE E at ok i te U aatgts sabi v 
‘ vad 1B bag ad stese 1 he Bisfee i Fe pel ete ho ata “ve 
Me een wicae a sagen My ASS ee vis ae - 
Lee S agbencee Sat * dadvicsd te Sk [Serra iitelits ‘ ‘ 29 we fil 
aaa Le ota hy ; ‘ + esis 428 eee eeek ; f : ah toc Mr Pe 4: RY: # 
OL Sa ns Ra f Us aes ee gy ist cy «ate SS as Mae 
‘ I 






+ 
os seme 















ar ae a Ll Bee | 
4 year sd 
RAD a caneneny aghye © 


ie ine O 
eae Sie 
Fug? 4. We 
ain S58 ehiee- Ld | ‘tents 


°. 4, 
at ads ig fe nh « pared EF) 





pee eat Fo af ae gtd 
Ot ot bt ae 
$ eeS 8 Ore ‘ 
at Ki "set get 
s Payee tas As . 
aed ie se 4 
o>, tee : 


~ 
ter tis wh Nelle | 
ay" (no ah PR a oot 
Ohne” a eqs 
Mt M48 OAM oop tefa, 
penne so seash shy 4ig! ; 
stele aves o 
Heuas Peete 
fap a eee 
bye yee Oe es 


t 
e 
af 


























PU oak oer ' bead 

































evs ald 









































Tgp reba " 


rang bests 




































































a mS eee b SANS +x Ses Fatale Wey 
Vee | | iae gel eo Ayre Ven =p . nen ty on 
tela g 2 e™ yee : oy ab rAades oosee Bet fee 
ere Poe rat 21% +n 10 detetee 8 oe (@ rptda Va gh ertet Seas . ore Lee at 
, Pe ona “hy wv ee \ Int | + gies “an ho kieia (*s moan a5 4 wipe 1 erent se ad pan 3 ru ‘ oN satis 
Alte bene Lak meee gngt at ov are ol Wie tee en yt ae Bek tale, ste fata” eabe 
te te ty iF Vette MGte tee Fe a1 ? ofa Fhe putas mab 5 Wace hea ‘tae satel ‘ ae Arr or i! avries: Ufa, rf : 
eu line CE UP IA eT sragset Abe ast mae eed ager t? Ase Sats (Syen fy, ait wing +, ea tetgs stis pa ey Maa RN 
Te Hela ALS Ne St ifs : SSRN henge pain OU ath het eone 8 aeyay alad gt oed peeled by vetet eat 
ie a irae toe "at tal Faned se TeN ee AL ogeet "tt ncbendseas aap eh tad eats! 6 A Liat teat ee MA aig) Gyths ita Weleda as 
nea atte ‘aNd ¥ ghee added wid bats ive Are his We dh, oA ME AALTICRL 4. af taboo shan ee Ab aE EN eth Th Wiz 
Me wee Fd Yoo 7 f esth ae | *f ats oy pia: th “ Mies 2M are ot) wb ls lef sfeeen) Se Ee sees 
aye etee Bi * a * a cor) Sa 
get a gee ade ieee : we hgh SA RRA GE a AG Me Sette tS 


eranctrigs aah yigae: aa, er errs 


bee perth 4 
7 wea La 
Mh Tosgea es 


pitha 
pt odd oy et 


ee? weapon PAS 





























































































































































































































































































































' “¢ ‘nt ete ie} es ¢ 3 2 Pf lene? 
Wh “gee Hy t tea Maly He heltet 1 ft z vad ite ‘ 
Ly we ahtat We bar 4 res ming Ae A ' 4 “4 t 4 we sv st 
=? ns s * , 
ia ie eae S 7 ht: 
ee | ee | ae GPE th ae! ae ete yet ee et ha te Sei We 
eee ar atdet © ditent ‘ ‘ ; cee ey vs tsi py a AS 5 MeAe 
Pig ‘geet ‘ ‘ Tle wet Tah adtae vod ror “ oe ee M3 Cea By Re me! anne Natit ae sf Rah tbe bY tin” Me tafe Pod tos en oe 
« «8 ‘ ‘ ° wees ee wee ph eae t 8G arg is ey peg ge fee te a*y f ougert paete ies ¢ 2. “ti Wea Ca deta | “gel %< a waft Se ers 
' ‘ Me Aa ee ak eer gene est Dar teat ie PECL mol Sake Seah Fat cate yee is : Ue Ese nae ? ace fra vay yeas a Sere Fenty ee ct she! Bote the fi 
. ie . * are t , Cy ead . id a Vt, 4 fen a a LS beet wgihye 4 waa” H « * s° r LL a + d ¢ a Rhea the oy Hy 
Ph as ; ieee fa. a os ‘ , <8 ” be SO t Hp ete Ae a . aals ia ‘ H tery ia ee ‘ 9 oer H yeteaee gt rgtet F114 Be gies tf" ah COS 3 les a. 
: ; ae a 4 Saker Bee eerie Puma, COME | Oe MC MI Ce Rk 7 ee ee ERP MEL AN qiat 
Ree an ane ’ Sr ireat eee alten ar wy ee Lee silat ete y at he oe vet AV egat bP Sisee 2 tame be >) en entesE Ng 
es once, ‘ ote ih Guerre foi wah tet VT ada iloe & 158, Vek Neh OOD aE hs ade O84 * isn a seats sh: 
ae . Ce eh I) t Soettise . . iz 1g ae ‘ cet hg: #3 AS aA 
, Hee : as 2 bys aes cis titeges ” 1 PER OGY ven xe shen Fe pe F iA 
ry ' cee H ! oof ‘ 4 “ete. 
aca hice fea Rau ee erat # Me baby hase tt yesiie te % ef aie es, “URARE 
Cane . sig ec Rip Gal ro aantg I Nes ts Saas bent ie gf San Pa) * sheng antsy F 
° : ") 8 “Spf etn me 18 Fa tee 8 Seth ‘ ; “ ve tt bs Shastri LE ak i etsy bed Ce v (Pur wigs. 
A Sree? rata yeerars AR Pere Yeas Wet ee we S de gegate f legs "thas WE oo Habe sdaa “eh ¢ ile, a Vee ~~! He oY 
et tue ' 1 er $6 mR ae (oe Vie actite @ of ee Ay pe ees EY wea dy 5 fs VED a ke De ws" 2 b¢y Jewett. 3 i A 
°% sone ; be oe Vere ye ag AU we ats tog EU | tOtd, AF at Ltebee he od! 4 oe hemes ter wee ar fo utiles aaa a ne 
to -ae so 1 . Aaa er ef Sha j ’ eo ts vie? Sr oe ee or 2 eT Haar site heey sy ete ky! #1 ol bet? odes 
’ F be . ee 6 A Ty ae ta hae Be eg te, aon ihe a Ca ches 5 fae Teh padad est ¢” eal a3 fn ryt aid 
ing Pi jo Fig h RIRS MEME ASNT RSE AMES GT 2 eae : Aare 4 a oak i vt siceelys hess “a fray . vd oe Be Bevis fee 
sete ti) Se td Vg ew We ds eee Me hee a ee NES e Re ane, fs +. esse beh ak, CAPS ota tard Seed 
. ; ’ bics: , re a ee “ POLY ee Sr . oo 4 a ‘s ee ‘ ‘ag Wty Mat * 4 Rien % Lar Oey Mole = fat 
WoeLe at Pay Ws ’ . Rae eee te lee oe th gt Dole cw) “he ne fie da ct as ohn Shes erg aut 7 os a Rete FL teary tits igs i 
r ray ' toy ae] Vives os : ah Hr fot btere a tC eats ba! ‘ 72 ts eee . 4 vy MNS ge iad Be wen 1 
A mo. sees ‘ ee bee 888 4 ferss Fe nae Bs tang Be a eevee ae IE Bid PA Ree iar theta vss 
‘ as alee H eenerg . white at pak ee jeg ry heidi atget! tah) PAE iar sir stg tleds sae 
+ eee ‘ Says > gis geo at 7 t “yee Sats aU USES: uy aes sca fas € 2 2 Py ae a cy. Pare) ot Wh sat Og ate tut, ae 
Doo4 . ee ace latte : yay : 1184 a ste a A ts qe arb ove Cire a ae ‘A Jeter se” tea toe Wake oa” aks hoe 
oi ap te bs Chat) tu iter UL Sai okie (iat ae, wil ” wie Mae) Slat ve® OF a “ate a? ‘; Py rae oe 
Cee fe a) ee. 8 re iy ae ere] austent Toe a) ate deista 4% AAAS: art atest ods ofata tye | 
a a | . . F Jae Yat haa | * . ' Pf gm PS wg hyhs Late? vemme ary: : 
aoe ’ . ’ ay tn ¥ ” u@ "Ray ke or fle Bie apes 
‘ 4 a ‘ oro Teer Parva ta | Ay cr ee } fess reo parker Te rey. ax 
toe ‘ ats Su Wceny etry erie PRE LA Mere aA ; is tens x 
bo tee aig oe feng Sh 4 sf Se Leatetee, eet os Did biota é as 2 Arata 2% ord 1 Whe 
«yw Pa ere a ery ie fers ob noe pee om t vipat ages se MEAL: wigs We Sy'y eiue Nay ale PEP Uat we 
a F Eat ae | vw! fe YF ose a satay? ot ue wea a > bebe 4! Mops heoeh be leas are A a *7 we 13 vice «. rh gt 
fut tay ' ‘ aa sty ieee t oe § ‘ “ a > ‘we, 5: PAV Fateh whe ‘ Veta te a eet ong f* Raat ‘ eye . $ oe * * Putas petst 
‘ ‘ . fers 2a lon Veae Pree tei Pre sa4 any aes seh ssa ft watt Neste she an Se ee “wer ise? otret 8 fs ay G! Vedwigatil ye « q at" “a stat tbare 
o 4 . oF 0 Gp oe * ont Oyta a fe Par ealess wren te RK OW Che oe fue tk y fey, feet “ . StePu te 048 aka gap Fstteg? qe eo gh 
en + re LQ tpt as st dag obs, ry a ie 2% rm 1° a5 sees wel om tty Ut ek Mee A A ed OS ae me cee ok Ly fe retest fice ete Sy eonse Resstgs dev, 
aaa | : se Fey ee) a eee a etaatae SCI PLAS GPO AC HE es Viates eb, ‘pli gtr gig Sow 3Fh 9 yeni 4% is 4 swan. weet Oo Cotas te ea 
vee * eo rs 44 te Ea pane ; Beane PY ee OL ' oe ob ' iva EO fa eh ee a CR re sytent at et 4 48, axe BPO PSrGr ts Boe Sage Ln Wie jet ned Chess 4 fae sgt 
- . pte : or ar : ty gt 5 te Gade a ew gg See 806 Se NaS Ue te! tinny . ae vagi tiene atte + ve Pe ct tet pra H vi daar ca rate Hee ite “" 4 es +) Brees iter ar 
eo ? ue . Ya ** fe ; oar oan q sh Be urty . . OE 4 « ae 
o) $ -ael a Trane oe fee gat we hd TE Wa ted t : Haste te en te ft ! vensetin ce 4 tae C Nos 3 phates 
tegot aek av ee SA \gheeud 4 
, es S Sire tall 2 
. weet , chia Ni doe t08 a 2. 
ae ye ener 3 ré et Naa, sar g ket oa * heyy? Apt Fs sine a 
Hees Mataayt "es Savy GS ih Let Mts, “Ws tis fie < vA wh o) 
ats te Cra vs ¥¢3 weed o.  Maye rm {3 ‘ 5 Th es Sah! wut ‘ 
rar eh VFS eagout as a Tag os, = gam itat a * ay ay tua te i Re A Ki 
1 (lame Pera Pe peers ates ghar. veedeln 2heees arse 2% £24! x Bele 
a hat reteset als pee wee sey Waet oy tty iain h2 ae iw 
eee. os ‘ uae er seghe, gly dee GN ee oe agin dé 
2,403 a 4 ’ a aS) $ Piny .s sista! 2 Sir) 7 sais 
} see we fur % $ IS soos wt Stig is 6 beans a Get 
Va *e 4.5 1M Persia x eee 
Or i e Mette ee ; Rea soi! 8 sae) ie 
tan C44 bt t , fa , ‘ Be Pte st 
‘ G as Es ‘ 
‘ ra | yhte 84 Hf ay 
nt " % 2 
stat oa Y 

























ae ae oy 


t b hd 
Reeehs 
it S0eiers 


‘ 
ae ee 
o. 
ee et tet te PO 
~ 





pai BE 34 





ary se 
ae 
- 
~ 
a en Ae 
Te ew wn an ky 
eh ft 








































































































































































a area et so hare) 
et 4 rrr a ' « ¥ 
Grecia: celat Re “ mee! ze agaete. BEATA x) tik ey Lae 
FD a | Cate te us, ee e @g.tee e8 hd” Lait Te ie “o Pedereee ow #1 ifade ste » ; 
aiden sritter dey : 4 Be wn pres Ber® ftesee® es ss : ‘ + 
sti dats tay ab at vasarg bee Gob haa Sees : serge Fea,® H 
vtaata gee 7 a 2 Me wiP det ohyt aye * Nb de itera + Pad th 1 “s. 
se fget Be Cokes Se se5lh tr Sei cnetyak fas Toe Brose aad atye? ‘ 
ion tana oe bets th Ss a Je 8 445 BOR ¢ 
BeOS de Tat tye gt ie amor G4terg® 2k FY " 
Woe Pe gk ase che Vaue® u s deayaye 4rd thee = z ay 
ore Saunt Wha at 4y fg 46 hate pred ste shew en: Rate got Right e Verge’ t 
: Fy Mot Cotes 7g fas ce oH wy Ht “aie ma as iene rhe PET aginst 
CY ay S., ur * : it rT . Pie oe, , 
Pete NEE tt amt Mo Set Mati a deeity eee Aiea ear porn ae 
ry ‘Lig artes * A a ac ne a ue uy yeit *f Piet ie. eer) tikte 's said 73 
PL %, if ed FSpeoreuaee ake oh ag ash ‘a ita se! Vika A He: yee tin ehh “a i* 
mie sats "Me a Cab be a Fd vee fy ernante tau* Ptarrase ’ Seal oye 
"4 gman atts 448 4 t mb aecesee gee Xgeey, ; 4S ts Ree* Stead of UT Fe eatee Ar tate ea 
#2) 6 see ett Nd ge its ri seg itaveteeta? 1 Pa oe Pek Pky a] 5} eer afer Fi Sait Si ats Stns yaa beh rips aay 
Beehee “7% De I ’ etged "oa t, 49! Py Sete ete eta PW Led ry eeu esas of state mc oetey a ‘es Gt rabie 
pate & Ww" , ° Fi . ar tale bic CRG is eta ts «Mette ft islet Sie oceagt Or a 
id eElaeey “haa ore Af eae nee Forge . ayy’ oe o 3 wate, Qo fee etereque y aehe Se Sak 
Dy ee ray ate? aot Sot pix tete Ode ot lag Pest oft bat 4% Vibe MES it wots ay lee™ here is 
“we ba eee ta! ’ qee asavat be toehes) eam be he Citar le Seger a be 
: Rous Tait Sasa Sites Sei sare Feagjenys 
a pRetee . rds os te dan eng wo Lg eS g° Uy 
WAL) We “gu ts At these? $y 445 Wa? an Ped he . 
» “ eos oot Cs Pid badass ¢ 4 
fat». tt « ’ APE erty heb Toh veh 
retey '¢. 6 tan @ id rat Ateahre ves “ 
WA OF Tensngesya } yo sy eeary, DEEN Arar 
“at Wa tes at ie d tas 54) Jas hae y 
eater ‘ 
$ | gk? oer ads : “eaeg al hy , 
jy vee eee ers. a . ‘ oases, 3 Pk ARMs : 
, 8 « ' wndy ht ge .Y) We pet letsse ? PeQeteebigue Baek, wid yigten 4, rhe , 
5; ; A he Sere £2 gre seeders gt * ffs eRSs TOF ee ds oY htc yer qed gtr Whe RLF bap 
ta > o teu F oe {we ate ’ PAs attr g oe bret neg .y iseeden, ea gtd > watt Mary » 
par vas sis" riper oe ewth ta eee Ap last owas 4 frei les rade rat 
‘oe ’ Ser palietta wt t a acs ee oe | »¥ aizd ee ynen wk au 
; a" A eee i Fh a toy rar, WAH #3 5° 5 Se 408 after ‘au ad 
. . . . HBr} ' Pa | ‘ . 
“ae . t. + 4, LP Fis tgs Ae vt jaan ce oh ane sia 1’t ‘yan y iF bed them t 
, v . fe Fe PLN aee : gery tt eer ate ap rats rbyead ee us 14e8at cg eae oS 
NS : ‘ is 7 ’ rr ‘ Canes eee ORI PCY ue 38 Pye eh es 
a 7 oo. oye a 4 “AKateee ty at 1 4A eceee bates e Val ts we ‘ ae i Delete iene 
' eon Tee ee ra ms his ' fw b onbe ge My yh floa oeee ae BTS cboee ee 4a *, ots + Tcckae, % 
aS . : sea antes : US ea Menace Siac tase Bee Teats tet feet % iv rie os 
. ts Datde af . a ” 
<i ’ wh ra 

























































































































































































































‘ Hy Lea |g joe ets pe tee? EN 
or ne Sere ih edi, states psete hitted by jae 
e "Vee oP Sp = 7 * ey 
ie Cag PT eee! CCE I 1 : Vy Ph oe “ 
; Ro guett aw the a Peg ts amet 8d at! : MS anes te ee "3 iets * 

aa eo est * OM ” 4 a teen egret og ate ve fitgewee mad? oat stages PR Ke sents it we 1% 

a? Bbstte out POP RT rt ogee € ite CT a abgheteh 928 gt! hee QUE Mop SUI, 444, ath € artes Prey) a 2 Ge ar em 

ety sg tyh oe AU ae send te ngaree gg wh ade ore (hbo! gk a taka AoE a tie eta 4 Dn ste 

bales ’ 63 aot oop ” % ia agcee no Pic) ' RP 3. Oi taie * + ties goes 
apy ele . 24 aes ngs the eu ttetpeaty & at ‘, Fas yr $j ww et eis 

: ae pre ’ gibes Re eee ‘ cas 
Ai © ee eee he 
“a Jor ate “ 2 egal a Sp Bgl ite Fhe ; 
, ie 5 Session El ny ays yd, ethos dated 

Bd Se mter . 4 
os nn ibe + ws of gases aie, 4 ' s ” read me - bit es 42 PTY a2 
ess . fa, ci lea oe? 1 rad) 4 oete Pe eG Phd 3. if sige iF) wey 
e .° s ; ‘ “ ery ms Aas > fy 

" "steely ere tetente ayes art ate wee NS , ts ce ton gig bel Fea i ay oe ‘f 
Z fer etnias BG Beit EB Maden a guard Ai ae Sih ae flee 
. 4. ‘ . ” te uJ t ae ial Se 5 : * = 
' ane ae : or rE rs oe Cee are ahi g Betas o yeades OLAS hae | ae’ 133° Ana's rr eres Wasi fake we ile: set fh et a edte a 8 
i ranks patie sets el are ee my pe ee a 18 cee ‘eats ae ue ret V4 Pd wart hy Messetse wai re puree Suatet 4S, seis aie a 1G aise ow 
Fy t . > beet Es ret sas = Z re tendin x u sae i 
owe ace 8 aE Pere ays ns gy ean sna PETS eee oes ye fea vet ie eae ¢ Fiyhe : 6 20s Ay ae AOE Bae grerayee 1 in laeehes the aH ule: aie { CP 
is raat ee “a ethoer Chama Segue megane ote Cr Cre eo ‘ ie UM Eis treet les eal Mee’ beta pean » OAT Erte oa it x tate *) ae Fak ah 
* = ts ote rag ees pide ete eee LI) ae fee rs Pe tae 383 wt ue PIR Fie faites won 2 ® at BG wa eth 8 6 ot oj Hor vy ee) Pee ei oye en 
praia ‘F ae wd " ed Bele tt as cen) S Ph Fd a LE ACE ICL oh Pt | ; vie eee “ st diss , fit ds Fe Ba pa Le artis | teat 
Ce al Hae tT Sa Ty Peat ee) : va Gat pet wisp e (ele te sleds sd Fe wd dm tee hie ; ie VK" {: pave v, a eynetts pha ty tir ¢ re 
‘ wt we « CBC ate og! as wt ef } er OO aie s od ' i) fs a ity oe ' wa Hiss, ra 1 Ky et ois seo an si gah a 
ae $ Paes dec Le eee eae rg ena Ly ay sar eyes os iy alas ‘it PA ty pana wrtates ies . * “¢ otha oka bet e, 
. ar ty ate Berk gs yay’ | ; ee aA War tenet ye eRe Ioege pty ee, 
; : : las T athe ee vat ad sipide Re aes 
i ae wpeboebe ey e23 i eens) Fe i he ty bose pay 
: ' oe iG qtisinenys ih. : nee eset 

Pa 3 antes laeateteatd Age (bee veule desta Sirs 
‘ ; tag te" 94 o04 neo ryt it ihges 

wpe id . 

“ ne ee “ uytigal are “4; wea if ae eae yt 
ete a fee Pek at, ation ai a) rage i TP Phot ear 
MG } ti tis Soe ae ra U ye ai th a ie oy tosancd Aste Rte ie A Cian bea Mb bent errs Lar 

. og fend tate Bua 8 + gs ‘9%, ae “eg OF a ne t eat ht etapenit Ae oe ak rid Ba hye 

' *7 ’ Pauses ’ ’ a6) LF eet) ef Ae eer | giety: ah ly ieee ang 3 ter fs wae i aE mf aie’ it 

: : c : = : - . ot Mater" trea Sys « be PEI Me Fgh ‘ty x io) 4g 5 . 
nels t: CA On SL “i nh titel sls $ ris . ret e of Bent pt . ih tat A iain 2 
7 ter het 


















a. ae Vege ty st 





Ree panel att 
ee saat gig ‘ Wn ay adh aah tiene a 
ue etg? Byars . eT heated 
‘ ne mt $5 oe vr 


eae 
ae < 


Pas bt Bae 
ena 

































erat 
ae resilais’ wt, reat 









































































































att Hae ne oss aetots 5 tice 
: se 4 Rion 
aE rei wely ie Cane % é artes nor 
. 7 y < 
; t ; nay w $e ee dats y tts irs nant . cent mess ee! 
YT op a Mato Bs eu Meike wala ete gy aye s qeardigey pictetee darter, ae 19 ho . what hos 
ee Tir ie aac berate bee iF pice aig Pelee tye PG a Petey Sclght nye alee iar 
TORN Winther 1K eta, sole fedh ute. 28 7 Rays pass ie gi sth mteig 
ek dee fan ' a yt. “i bs : Uaytyez, ag a 
DV oieret were Saye a Ce Dg f iad eS as 
Ue ed a : 4 intets bit ee 
: rt = "a fe gy at a a mye mt he see ae gigs 
4 yy" Sede s op er ey 
H r : auth sideaece Sue pyrene 
: at RESO RR E eS What eth shia et 
Pay oe Meow seiagiete grant Maines oy 
SUL Aroe ALS i ie Rt oes, 
patterns 
vat ae «" A fase? 





pee S| Yaw a 





‘ «fey eins 
ut vie itd Pe me Anes 
ROA aN eR NTE fe 
; A wy is wp " a vied ite eiaae be 
nau a veh tite ds do! “ah UB ened eee ee 
Pip Or ww ry a 
<j. oy) ter ines, & xs Ws eh ee Ye ty iets ” 

SAT ace! 
2 ats fram 



































ae er 
sate rts 
ie et Re 


ie seo ay 





as oe Dae, * we 
Yabed mae n 
av 85 stent "he 
4 fgtg te! 1 oy 
+t wl d 



















: ' 4 Te tess30 cle 
1 ' whe re . she ded 

oan $i . rai UR i Sesh es, oan br bh di ataa ] 
98> 5 Na eRe V1 tats 6 lye S padey BakG. w apiste 


DUDLEY KNOX LIBRARY 
NAVAL POSTGRADUATE SCHOOL 
MONTEREY CA 93943-5101. 











Approved for public release: DiStribution is unlimited 


Orbital Maintenance of Endoatmospheric 
Low Barth-Orbiting Satellites 


by 
paved) D. Pauls 


Lieutenant, United States Navy 
B.A., Physics, Eastern Washington University, 1982 


Submiuttecein Pareeral foalfiliment of the 
requirements for the degree of 


MA ener 6 oC il ENCE 
IN ASTRONAUTICAL ENGINEERING 


from the 
NAVAL POSTGRADUATE SCHOOL 


DECEMBER 1991 


Unclassified 
SURITY CLASSIFICATION OF THIS PAGE 


Form Approved 
REPORT DOCUMENTATION PAGE OMB No 0704-0188 
1. REPORT SECURITY CLASSIFICATION 1b. RESTRICTIVE MARKINGS 
Unclassified 
1. SECURITY CLASSIFICATION AUTHORITY 3. DISTRIBUTION/AVAILABILITY OF REPORT 
Approved fox publike release: 


FICAT . ; os: 
.. DECLASSIFICATION DOWNGRADING SCHEDULE Distribution is unlimited 








PERFORMING ORGANIZATION REPORT NUMBER(S) 5. MONITORING ORGANIZATION REPORT NUMBER(S) 
1. NAME OF PERFORMING ORGANIZATION 6b. OFFICE SYMBOL 7a. NAME OF MONITORING ORGANIZATION 
(If applicable) 
VYaval Postgraduate School AA Naval Postgraduate School 
;, ADDRESS (City, State and ZIP Code) 7b. ADDRESS (City, State, and ZIP Code) 
Monterey, CA 93943-5000 Monterey, CA 93943-5000 
|. NAME OF FUNDING/SPONSORING Bb. OFFICE SYMBOL 9 PROCUREMENT INSTRUMENT IDENTIFICATION NUMBER 
ORGANIZATION (if applicable) 
ADDRESS (City, State, and ZIP Code) 10. SOURCE OF FUNDING NUMBER 
PROGRAM PROJECT TAS WORK UNIT 
ELEMENT NO. NO. NO. ACCESSION NO 





. TITLE (inciude Secunty Classification) 


ORBITAL MAINTENANCE OF ENDOATMOSPHERIC LOW EARTH-ORBITING SATEDDEUiS 


. PERSONAL AUTHORS 
AVEID Dee Ais 


a. TYPE OF REPORT 13b. TIME COVERED 14. DATE OF REPORT (Year, Month, Day) 15. PAGE COUNT 
laster) sie sme PSO Oa DECEMBER 1991 74 
. SUPPLEMENTARY NOTATION 


he views expressed are those of the author and do not reflect the official policy 
Xr position of the Department of Defense or the U.S. Government 


. COSAT! CODES 18. SUBJECT TERMS (Continue on reverse if necessary and identify by block numbers) 
1D GROUP SUB-GROUP orbits, satellites, | low-Earth orbiting 
ses a SS satellites, fuel-minimization 


. ABSTRACT (Continue on reverse if necessary and identfy by block numbers) 

The optimization of spacecraft trajectories in vacuum has received extensive consideration 
ince the inception of space flight, yet, the effects of atmosphere have been largely neglected. The 
dvent of low Earth-orbiting, large satellites and platforms necessitates that atmosphere be included 
n the optimization process. A practical means of studying this topic is as a problem in minimum—fue@ 
rbital maintenance. Optimal control theory advances the notion that orbital maintenance is optimized 
hrough periodic tnrusting as opposed to forcing Keplarian motion by nullifying the effects of drag 
th thrust. Furtner, this must be optimized by primer vectoring. This thesis examined the efficienmm 

a simple method of orbital maintenance using fixed-angle transverse thrusting. Results show that 
or the purpose of fuel-minimization, the width of tne radial band in which the satellite is toume 
aintained, is dependent upon thruster size. In nearly all cases, a thrust-angle of 70 degrees 
aximized the fuel saved. Tnis thesis shows that fixed-angle transverse thrusting does not improve 
n forced Keplarian motion and hence thrust vectOring mUSt Seucpein zea. 


J 
c 
— 
‘Ss 
-_ 


DISTRIBUTION/AVAILABILITY OF ABSTRACT D1. ABSTRACT SECURITY CLASSIFICATION 
XX UNCLASSIFIED/UNLIMITED __ SAME AS RPT __ DTIC USERS unclassified 
3. NAME OF RESPONSIBLE INDIVIDUAL P25. TELEPHONE (include Area Code) | 22c. OFFICE SYMBOL 
I.M. Ross (408) 646-2716 AA/RO 
Form 1473, JUN 86 Previous editions are obsolete. SECURITY CLASSIFICATION OF THIS PAGE 


S/N 0102-LEF-0Oi4—Gcts Unclassified 
aL 


NAVAL POSTGRADUATE SCHOOL 
Monterey, California 





ORBITAL MAINTENANCE OF ENDOATMOSPHERIC 
MOWSEARTH-ORBITING SATELLITES 


by 


David D. Pauls 


DEGEMEpE Ro 1. 


Thesis Advisor: 





Approved for public release: Distribution is unlimited 





ABSTRACT 


The optimization of spacecraft trajectories in vacuum has 
received extensive consideration since the inception of space 
flight, yet, the effects of atmosphere have been largely 
meqgrected. The advent of low Earth-orbiting, large satellites 
and platforms necessitates that atmosphere be included in the 
optimization process. A practical means of studying this topic 
is as a problem in minimum-fuel orbital maintenance. Optimal 
control theory advances the notion that orbital maintenance is 
Sepeimized through periodic thrusting as opposed to forcing 
Keplarian motion by nullifying the effects of drag with 
mij~iict. Further, this must be optimized by primer vectoring. 
This thesis examined the efficiency of a simple method of 
Orbital maintenance using fixed-angle transverse thrusting. 
Results show that for the purpose of fuel-minimization, the 
width of the radial band in which the satellite is to be 
maintained, is dependent upon thruster size. In nearly all 
cases, a thrust-angle of 70 degrees maximized the fuel saved. 
This thesis shows that fixed-angle transverse thrusting does 
not improve on forced Keplarian motion and hence thrust 


vectoring must be optimized. 
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NOMENCLATURE 


radial acceleration 

transverse acceleration 

semi-major axis of an ellipse 

negative reciprocal of thruster exhaust velocity 
Seeetricient of drag 

magnatudesof drag 

radial component of drag 

mransverse component of drag 

total energy 

Seecernal forces in radial direction 
external forces in transverse direction 
Gravitational acceleration at sea level 
Pom@emyvagin’s H-function 

specific impulse 

Mere, VECeor in Girection of thrust 
spacecraft mass 

mesition vector 

radius 

maximum radius 

minimum radius 

initial radius 

Zmomal radius 

path travelled by spacecraft 

iii¢tial position of spacecraft 

final position of spacecraft 
aerodynamic reference surface of vehicle 


ilmrust vector 

magnitude of thrust 

Gadial component of thrust 
transverse component of thrust 
time 

initial time 
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velocity 

VeElOCi Ey Wee tors 

radial component of velocity 

transverse component of velocity 

thrust angle 

exponential density scale factor 

change in total energy 

radius at which periodic thrusting is commenced 
change in time 

specific energy 

orbital coordinate used to define spacecraft position 
Lagrange multiplier 

Mess .COsta ce 

velocity CoOSteteu, ere 

radius costate vector 

Gravitational concwane 

atmospheric density 

atmospheric reference density 
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I. INTRODUCTION 


mm. 1963, Dees ibawden, i” Mls m—nemocraph [|Ret. 1], 
laid the foundation for what has become a highly sophisticated 
S@Ooiscipline of astrodynamics, optimization of space 
Meaeyectories. In 1979, Marec [Ref. 2} provided a more 
comprehensive treatment of the subject in his text on optimal 
meremeenm trajectories. Examination of the optimization of 
Spacecraft trajectories has been treated by many authors in 
manners similar to these two great works, yet, until recently, 
the study of atmospheric effects upon these trajectories was 
largely neglected. Research and development of hypervelocity 
vehicles have kindled interest in this area, through which 
Study, other areas of interest have emerged. 

One such area is the effect of aerodynamic force on non- 
Premio, Or, Dlunt bodies. First addressed by Ross and Melton 
(ker. 3), this subject is of particular interest for 
two reasons. First, atmospheric effects on low-Earth orbiting 
(LEO) satellites are of obvious interest and second, as stated 
Mmeedeir paper [Ref. 3:p. 2), better understanding of this 
phenomenon could provide deeper insight into the more 
Semeiicaced topic of lifting bodies in the upper atmosphere. 

Much work has been accomplished in this area pertaining to 
Meemrare prediction of satellite orbits (kee. 4, 5)" 


The focus on atmospheric effects as they pertain to the 


specific problem of minimum-fuel orbital transfers, however, 
1S UNIQUE temic «| and the follow-on work described here. 
The particular problem of orbital maintenance can provide 
insight into the more general, and complex problem of orbital 
transfer. 

Historically, orbital transfers (coplanar) have been 
accomplished as either, two, or, three-impulse maneuvers 
[Ref. 6:pp 78-88]. The problem of orbital transferee 
approached as a minimization of energy required to move a 
satellite from orbit A to orbit B, or, equivalently, a 
minimization of the characteristic velocity. A satellite that 
has descended from an initial orbit due to a disturbing force 
such as drag, and which must be returned to its initial orbit 
can be approached as a problem requiring an orbital transfer. 
Orbital transfers such as this are optimally performed by two- 
impulse transfers, such as a Hohman transfer. Large orbital 
transfers (eee are optimized with three-impulse 
maneuvers [Ref. 6:p 87). As low-Earth orbits become more 
frequently utilized, deeper understanding of the effects of 
drag must be achieved in order that orbits, propulsion systems 
and costs are optimized. This is particularly applicable for 
large satellites, such as the proposed space station, that 
must remain in low orbits for extended periods of time. 

Additionally, more complex areas of study, such as that of 


lifting bodies in the upper atmosphere, could benefit from the 


insights gained through a deeper understanding of atmospheric 
Serectes ON Non-lifting bodies. 

Research on the atmospheric effects on low-Earth orbiting 
Spacecraft is Sure to receive much attention in the future. 
The benefits to existing and future systems, while extensive, 


remain relatively unexplored and demand the attention of the 


aerospace industry. 


A i GENERAL FORMULATION OF THE PROBLEM 


Optimization of orbital transfers is a subject thai 
achieved a high degree of sophistication and many elegant 
solutions exist [Ref. 1, 2]. However, the specific treatment 
of non-lifting bodies is in the initial stages of development. 
The increasing number of low Earth orbiting satellites 
requires that a study of atmospheric effects on orbital 
trajectories be conducted. In this thesis, the problem of 
minimum fuel orbital maintenance is considered. Two methods 
are examined by which orbital maintenance may be performed. 

One method is to counter drag with thrust. In this fone 
Keplarian motion, the reaction control system would thrust 
continuously for the duration of the satellite’s lifetime with 
magnitude and direction equal and opposite to drag. The second 
scheme considered here, utilizes periodic transverse 
EHrUSE ne, Oe, non-Keplarian MOELlom to cerrec: for 
perturbations due to drag. The question to be answered is 
whether or not optimal non-Keplarian trajectories are superior 
to forced-Keplarian trajectories. 

Let the problem a defined as maintaining an orbital 
deviation within a specified radial band r,,,SrSr,,,- Is forced- 
Keplarian orbital maintenance, i.e., exactly countering 
aerodynamic forces with thrust, superior to non-Keplarian 


orbital maintenance, i.e., allowing the orbit to decay to a 


certain point, then reboosting to a point above the desired 
eeertuce. While maintaining the spacecraft within the 
specified radial band? The study performed by Ross and Melton 
[Ref. 3:p. 4] suggests that forced-Keplarian motion is not 
optimal and that thrust vectoring must be considered if an 
optimal solution is to be obtained. 

This thesis addresses an additional question, can orbital 
maintenance be optimized if thrusters are fired at a fixed 
o@epee tO the local horizon and if so, what is the angle or, 
preferably, range of angles that achieve optimality? 

Ross and Melton [Ref. 3] develop their theory through the 
methods of optimal control theory. It is proposed here that, 
while ideally accomplished in that manner, satisfactory and 
enlightening results may be obtained through the use of 
relatively unsophisticated mathematics and the aid of computer 
modelling techniques. To exemplify this statement, let us 
examine the mathematical development in Ross [Ref. 3: pp. 


ieee Drag iS given by 
Pare eg | 
D= “ape CO ae (1) 


where p(r) is atmospheric density, S,.. is the reference 
Surface area of the spacecraft upon which the aerodynamic 
Mem@ees act, C, 1S spacecraft’s coefficient of drag and v is 


Meme tocity. The equations of motion can be written as 


n= Vv (2) 


eee ———— py, Trex) (3) 
r?3 2 m m 
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Applying the principles of Pontryagin, the Hamiltonians 
given by 
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The costates are then written as 


es = OH = Bieter oe v\ | 9p) _ aes + = 6 
a, = HE = {Ssec0o) av (4) (32) - Ha sua (S)e © 
Pav (2) (7) 
m V 


_ 2 OH = Sret“p al 2 oe 8 
An = ar | aal pve (+) Tradl (=| (8) 


A closed form solution to these equations does not exist and 
only after initial and boundary conditions have been 
determined may the solution be obtained through numerical 
integration. This is a cumbersome method for determining the 
optimal direction of thruster firing. In this thesis, woemiegs 
into the possibility of a constant vectoring scheme that may 


result in nearly identical performance. 


A computer generated model will strive to maintain 
spacecraft trajectory within +1.5 kilometers of the injection 
altitude utilizing a periodic fixed angle transverse thrusting 
control scheme whose direction iS maintained at a fixed angle 
relative to the local horizon. The propellant consumed will 
then be compared to that of the same satellite employing a 
control law that sets thrust equal to drag at every point 


[renin the orbit. 


III. SIMPLIFIED PARAMETRIC FORMULATION 


A method less elegant than optimal control theory, but 
nonetheless valid, is that of parametric variation. The 
equations of motion are developed for the satellite’s orbital 
trajectory and certain parameters varied to achieve "optimal" 


control of orbital vVariae tens. 


A. DEVELOPMENT OF THE EQUATIONS OF MOTION 

For the purpose of building a solid foundation, Gemeaums 
Simplifying assumptions are made. Orbital motion is assumed to 
be coplanar, the initial spacecraft orbit is aSsSsumeqdiaa 
circular and since the spacecraft is a non-lifting (Sig 
body, drag is the net aerodynamic force acting Upemeaee 

The equations of motion for this two-body system can be 


written as 





ae 
a, — Dy ar (9) 
aoe a0 


where a, is the radial acceleration of the spacecraft, a,, is 
its transverse acceleration, F, is the sum of the external 
forces in the radial direction, F,, is the sum of the Jexteiee 


forces in the transverse direction and m is spacecraft mass. 





The external forces acting on the spacecraft are the 
Gravitational field, aerodynamic forces and thrust. Figure 1 


illustrates the coordinate system and the net forces. 


Transverse axis 


Satellite orbit 


Figure 1 Graphical Representation of Coordinate System 


Referring to Figure 1, it is clearly seen that the components 


for drag are given by 


S 
i" 


-D sin(y) (11) 


S 
" 


-D cos (y) (12) 


Likewise, thrust is written as 


T, = Tae Oouay (73) 
T., = 13& Stray (14) 


The angle y is the flight path angle, defined by tie 
intersection of the velocity vector and the transverse axis. 
The angle @ is called the thrust angle, defined by the 
intersection of the thrust vector and the transverse axis. The 


equations of motion, then, become 


2 eae (15) 
re m m 


GO; + 207 = =-22...— (16) 








where HW is the Earth’s gravitational constant. 


B. DETERMINATION OF THE CONTROL VARIABLE 

Following the development of the previous chapter, it is 
desirable to maintain the spacecraft within a radial band of 
a predetermined width. This may be accomplished by directly 
controlling either radius or eccentricity. Radius "ise 
control variable of choice for a number of reasons: jigs 
found directly from integration of the equations of motion, 
changes are easily visualized and radius control provides 
indirect control of the eccentricity. It is clear Cnaniias 
varying the thrust, control of satellite radius is possible. 


Examination of the thrust equations (Equations 13 and 14), 


AL 18, 


presents two methods by which thrust may vary, changes of 
amplitude or change in the direction of the thrust vector. 
The following chapter presents a method by which radial 
@Mewration is controlled through variation of the thrust angle 
and then tested by varying the thrust magnitude. A computer 
model is developed that simulates the trajectory of a 
Spacecraft, graphically demonstrating the effects imposed on 


@emeemrough variation of the direction of the thrust. 


ude 


IV. DEVELOPMENT AND TESTING OF THE COMPUTER MODEL 


A. COMPUTER PROGRAM DEVELOPMENT 

As indicated in the previous chapter, a computer program 
is developed to simulate spacecraft orbital trajectories and 
1s included in Appendix A. The program is written in the 
fortran language and employs a fourth-order Runge-Kutta 
numerical integration routine to integrate the equations of 
motion. The program consists of six sections, the main program 
and five subroutines. The main program controls input and 
output while the subroutines provide various other functions 
necessary for accurate simulation of orbital trajectories. 

The first subroutine calculates drag experienced by the 
spacecraft. Ijatstam i= liail 7 a model incorporating constant 
atmospheric density is used. This facilitates verification of 
the program, after which, an exponential density model is 
used. It is acknowledged that more accurate atmospheric 
density models exist, however, the exponential model provides 
Satisfactory accuracy over the range of travel experienced by 
the satellite (+ 1.5 kilometers from initial orbit R,) as 
controlled by the simulation. The second subroutine contains 
the equations of motion governing the spacecraft’s orbital 
trajectory. To facilitate handling, the equations were broken 


into parts. Solving Equations (15) and (16) for acceleration, 


ee 


1t is seen that the right-hand side of the radial equation has 
four components and the angular equation has three. These 
components are labeled A, B, C, and E for the radial equation 
pee) and S for the equation governing angular motion. The 
@mmree Subroutine contains the fourth-order Runge-Kutta 
numerical integration routine used to integrate the equations 
of motion. The next subroutine calculates the parameters of 
the satellite’s osculating orbit. The last subroutine contains 
the control law governing the activation and deactivation of 
the thrusters responsible for the periodic maintenance of the 
Eeaieellite’s orbit. 

The following parameters define the specifications around 


which the computer model was developed. 


* Spacecraft mass = 20,000 kg. 
@eepacecraft frontal area (S,.-) = 60 m*. 
meeoertficient of drag (C,) = 2.2. 


Meeititude of Spacecraft’s orbit (R,) = 260 km. 


* Atmospheric density at R, (p) = 8.3130 x 10°" kg/m’. 


B. PROGRAM VALIDATION 

Program development proceeded in stages, with each stage 
Peemaring validation prior to beginning the next. Initially, 
all external forces except gravitation were neglected. 
Clearly, radius, speed, angular momentum and specific energy 


must remain constant for the program to be considered valid. 


1S 


Having accomplished that elementary stage, phase two 
introduced drag (Equation (1)). Taking into account themieaee 
that aerodynamic forces acting on a spacecraft are very small, 
and consequently, the time required for significant changes to 
occur, very large, atmospheric density was increased by 
approximately three orders of magnitude in order to reduce 
computer run time. 

Initially verifying that radius continually decreased, the 
accuracy was tested by comparing the difference in altitude 
per orbit calculated by the program to that Caveumeiee 
manually by simplified equations. This is accomplished through 
the use of a "rule-of-thumb" and is developed in the following 
section. 

1. Development of a Rule of Thumb 

Work done by drag is a function of the path travelled 
by the spacecraft. Therefore, the amount of work done 
corresponds to the change in energy of the spacecraft, which 


1S given by 


Work done = AE = Drag x 2nr (17) 


Spacecraft specific energy iS given by 


g- vi _ be (18) 
2 rg 2a 


Total energy E is equal to the specific energy multiplied by 
spacecraft mass. Assuming that the spacecraft is in a circular 


orbit, the semi-major axis a is equal to the radiusiime 


14 


Performing this substitution then differentiating both sides 
of the equation while holding mass constant yields 
de = #@ ar (19) 
Darr 
meme Fouations (17) and (19) equal to each other and 
solving for the change in radius yields 
ia ge (20) 
Lm 
Equation (20) is the decrease in radius per orbit of the 
spacecraft due to drag. Despite inaccuracies resulting from 
the simplifying assumptions, this rule of thumb is accurate to 
within a few percent. The difference in spacecraft radius 
calculated by the computer program matched that of the thumb 
rule within a few percent thereby validating the model through 
this point in its development. Phase three introduced thrust 
while setting drag equal to zero. Clearly, any results other 
meee StCadily increasing radius, angular momentum, and 
specific energy would have been cause for program 
invalidation. 
2. Development of a Control Law 

The purpose of the control law is to maintain 
satellite radius within a prespecified bandwidth. By 
monitoring certain variables, activation and deactivation of 
Station keeping thrusters can be determined. Keeping the 


control law simple in order that computer memory and run time 


Ve) 


related to station keeping be kept to a minimum is alsoma 
desirable goal. With these facts, control law design proceeded 
as described in the next section. 
3. Control Law design 
a. Approximate method 

Specific energy is a function of spacecraft Pages 
For that reason, specific energy is a very useful parameter 
that can be used to maintain the satellite within the 
specified bandwidth. The energy lost when the satellite’s 
radius decreases, is dependent upon the path taken by the 
Satellite in its descent as illustrated in Figure 2 below 


(exaggerated for clarity). 





Figure 2 Path Travelled By A Spacecraft Between Two Orbits 


If S represents the path travelled by the spacecraft, then 5S, 
is its initial position and S, its final position. Thee aiaee 


acting against spacecraft motion is drag, which directly 


16 


Meeeses the velocity vector. With this in mind, the energy 


loss iS given by 


S2 


AE = [> ds (21) 
Sy 


Miemarc length S is a function of radius and angle turned 


moeeeugi, and is given by 


Ss = 76 (22) 


Differentiating the above equation yields 


adS = rd0 + 6dr (23) 


P@leemrtuting EQuation (23) into Equation (21) yields energy 


Toss in terms of the known variables, r and @ 


Rr 6 


= (24) 
AE [aids + pre 


@memeying Simpson’s rule, the loss in energy can be 
approximated fairly closely. Assuming atmospheric density, 
flight path angle, thrust angle and thrust are all constant, 
Simpson’s rule applied over ten iterations yields the energy 
eeeewin moving from R, to R; 


ae =| 
30 


a2 (25) 
|pe B80) + 30 PE (29) 


femerting tO Figure 2, R,; is the point at which the thrusters 
feeerrre. In terms of the control law, R; is the control 


variable and can be relabelled as AR Since it is variable and 


hey 


determines the width of band in which the satel iaigeuaee 
maintained. As developed for this model, AR is the initial 
radius R,; minus one kilometer. This value is chosen so as to 
maintain maximum orbital deviation within + 1.5 kilometers of 
Ro. Now that the energy loss has been approximated, the 
objective is to determine the length of time to fire the 
thrusters in order to replace the energy. The return path of 


the spacecraft is a function of the thrust. 
Ag = {(T-D) + as (26) 


Recognizing that ds is related to the time rate of change of 


the position vector or arc length S, Equation (26) becomes 
Ag = [(f- Dvdr (27) 


Reducing Equation (27) t© Competent weerm ly te lce 


Ge ce 
AF = | (7,-D)iv.dt | (i eae wea (28) 
fo fo 


Integrating and solving for At yields the length of time that 
the thrusters must fire in order to replace the energy iost 


due eordrac: 


At= aE 


a 29 
ne 7 D,) Vv, + (ie ~ D.,) Ver oi 


Examination of Figure 1 reveals that velocity can be written 


as 


18 


Ve = VST Eye (30) 


and 


V,, = veoos(y) (31) 


r 


Potices Ieee nrough SGewer  AppendsceeBe illustrate 
results obtained from this method for thrust angles of 60° and 
Peeeeexamination of plots of thrust versus orbit (Figures 2 
ana 5) reveal the inadequacy of this model. Thruster firing 
times are seen to be on the order of orbits rather than 
Met@eons Of orbits. This is due to two factors, first and 
most important being that the energy change is calculated from 
the initial orbit to the point AR where the thrusters begin 
firing. The problem arises from the fact that the satellite is 
Still in a descent at this point and continues to descend 
until its motion is reversed through the opposing force of the 
thrusters. As a result, the satellite loses more energy than 
is replaced. The second problem arises from the inaccuracies 
inherent in the assumptions required to perform the 
approximation. While the first problem is easily resolved, the 
changes in atmospheric density and flight path angle, while 
very small, are not constant and the resulting errors combine 
Pemeciaer this model unsatisfactory. A more accurate method is 
to calculate the energy loss directly using the variables 


ea eived #rom integrating the equations of motion. 


ee, 


b. Direct method 

Using Equation (18), the initial and instantaneous 
energies may be calculated at any instant during the 
trajectory of the spacecraft. As before, thrusters will fire 
when the spacecraft descends below AR. The program then 
calculates the spacecraft’s specific energy each time the 
equations of motion are integrated. Comparing this value to 
that of the initial circular orbit, the control law Gomme 
the thrusters to fire until they are equal. Examination of 
Figures 1 through 6 of Appendix C reveals that while results 
are closer to those expected, this method also appears to be 
inadequate. Thruster firing times (Figures 2 and 5) reieus 
excessively long. Reevaluation of the model suggested that the 
solution might be a function of the thrust to drag Yama 
IncreaSing thrust to a value of 300 N and then Cxaminmnag 
results for thrust angles of 60°, 65° and 70° revealed this to 
be the solution. Results are found in Appendix D. Thrust 
angles of 65° and 70° maintain radial deviation withammieae 
prescribed bandwidth of three kilometers with a trend that 
indicates they will remain so indefinitely. Thrust plots 
(Figures 2, 5 and 8) illustrate that the thrusters are firing 
over a small portion of an orbit instead over a period of many 
orbits as before. Plots of spacecraft radius versus OrbiUgeme 
included in each appendix to illustrate the success of the 
control law in maintaining radial deviation within the three 


kilometer band. Energy plots are also included as 
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Semeeoberating illustrations of the Spacecraft’s energy level 


at €ach point in its trajectory. 
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V. ANALYSIS AND RESULTS 


Analysis of the validated model is performed in four 
steps. Prior to testing, drag was returned to its true value 
of 8.3130 x 10°** kg/m’? and thrust was set equal to a valuemum 
25 N. This value is determined by multiplying the thas aia 
drag value of the test model by the drag actually experienced 
by the satellite. The thrust to drag ratio with the increased 


densaty is 


300N 
——— = 75.834 32 
3. 9SGry (°a) 


With 9 = 8.313 x 10°** kg/m’, drag is calculated to be OMogmmr 
Multiplying this to the value in Equation (32) yields a tiga. 
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A. USE OF A CONSTANT DENSITY MODEL 

Initial testing maintained a constant density while 
varying thrust angle. The results obtained from these tests 
are interesting. As illustrated in Appendix E, approximately 
Six and a half orbits are required for the spacecraft’s radius 
to decay to AR. Thruster firing places the spacecraft into a 
Slightly elliptical orbit (typically, eccentricities were 
found to be on the order of 107°). Examination of energy plots 
reveal that the energy of the elliptical orbit is very close 


to that of the initial circular orbit. This results wee 
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large energy change as the spacecraft initially descends from 
Ro and is returned, then, subsequent small changes are 
required once the spacecraft is established in its 
Seweyptical” orbit. A thrust angle of 70° is seen to yield the 
best results (see Figure 6 of Appendix E). Radial deviation 
remains just within the specified radial band with a trend 
that indicates it will do so indefinitely. At angles of less 
than 70°, radial deviation steadily increases until it exceeds 
the prescribed limits. An illustration is provided in Appendix 
E, for a thrust angle of 65°. Thrust angles greater than 70° 
eeumeeow a trend illustrated by Figure 7 of Appendix E (a@=75°) 
until reaching approximately 85°. Above 85°, increasingly 
larger values of thrust are required to maintain the 


spacecraft within the prespecified radial band. 


B. USE OF AN EXPONENTIAL ATMOSPHERIC MODEL 

Use of a constant atmospheric density model permitted the 
determination of an optimum angle that could be compared with 
that determined by a more accurate exponential atmospheric 
density model. Also, realizing that density changes would be 
small within a three kilometer band, the constant density 
model divulged a reasonable facsimile of results obtained from 
the exponential model. As stated previously, more accurate 
density models exist, for example, the Jacchia Atmospheric 


Model (J70). The exponential model, however, provides 


a 


sufficient accuracy for the development and testing done in 


this thesis. With this in mind, density iS now Givermmee 


p = p,e bir) (33) 


where Bf is determined to be -2.12 x 10° m”’ [Ref. 7], 
As predicted, the resultS are nearly identical to those 
obtained from the constant density model. A thrust angle of 
70° maintains the spacecraft within the three kilometer band 
with the same trends as described in the previous section. 
Results are illustrated in Appendix F. 
1. Comparison with a Forced-Keplarian Model 

It is now possible to compare results obtained from 
this model with that of a spacecraft experiencing forced- 
Keplarian motion (thrust equal to drag, resulting in a 
Circular motion). Modifications to the program t0O Olmeansauee 
model in which thrust is equal to drag are very simple. The 
drag subroutine calculates drag then sets thrust equal stoma 
The subroutine containing the thrust control law is removed 
from the program entirely since thrusters will fire 
continuously as long as propellant is available. Drag always 
opposes the velocity vector which, as previously shown, is 
defined by the flight path angle yY. If the thrust 1S xen 
opposite and equal to the drag force, then the thrust angle @ 
is equal to the flight path angle. Therefore, modifications 


consisted of setting @ = y, removing the subroutine containing 
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mene nrust control lawwand setting”™thrust equal to the drag 
Calculated in the appropriate subroutine. Since the spacecraft 
Meera cially in a circular orbit, validation of these 
modifications is achieved when the spacecraft’s orbital 
parameters remain unchanged over the test period; in this 
eaeen 20 orbits. 

Finding that the model performed as expected, plots of 
propellent consumed over the test period are compared with 
those for the spacecraft experiencing non-Keplarian motion at 
a thrust angle of 70°. Results are contained in Appendix G. It 
is clearly seen that orbital maintenance using forced- 
Keplarian motion is superior to that uSing non-Keplarian 
Monon . 

ee Further Testing 

To determine the "robustness" of the model, two 
additional tests were performed. In the first, thrust was 
maintained at 25N while specific impulse was varied over the 
range of 200 sec, valid for chemical reaction engines, to 2000 
sec which is valid for electrothermal engines. In the second 
test, specific impulse is returned to its initial value of 300 
sec while thrust is varied over a range extending from 1N to 
Sone 

a. Constant Thrust, Varied Specific Impulse 

Varying specific impulse while holding thrust 


Pemereant, forces the rate of fuel consumption to change. 


eo 


Specific impulse is given by given by 


J oy. Sle (34) 


where I,, is the specific impulse, T is the magnitudem= us 
thrust, g, is the gravitational acceleration at sea-level and 
mis the change in mass. 

Clearly, decreasing thrust results in a corresponding 
decrease in the rate of propellant consumed over a given 
period of time. This is graphically represented in Appendix G. 
The absolute quantity of fuel consumed decreases with 
increaSing specific impulse, but the percent difference 
between the reboost and forced Keplarian methods remain 
virtually unchanged. 

The percent difference is calculated by taking the ratio 
of the values of propellant mass at a specific time for 
thrust-equal-drag plots and reboost plots to determine the 
relative change between them. This provides a truer comparison 
of the two schemes than does simply comparing the end values 
of the plots. 

The results illustrated in Appendix G indicate that, while 
playing a Significant role in propulsion system optimization, 
specific impulse is not a function of the method used to 
perform orbital maintenance and will not affect the particular 


outcome of one method more than another. 
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b. Constant Specific Impulse, Varied Thrust 

The previous section illustrated the significance 
played by specific impulse in minimizing fuel expenditure 
during orbital maneuvers. Obviously, engine "Size" plays an 
equally important role in that process. It is expected that as 
thruster size increases, fuel expenditure will increase. 
Appendix H bears this out. While it is theoretically a simple 
matter to choose the proper specific impulse (bigger is 
better), this is not the case when choosing thruster s1Ze. 
Examination of the plots in Appendix H reveal that when 
comparing engines over a certain time span, the end results do 
Memmeprovide a ready solution. Figure 1 is a case in point. 
meee thas case (1N thruster, I1,,=300 sec) results in the 
least amount of fuel expended over the range of thrusts 
@iiesen, it is obvious that this is not a wise choice of 
Mmewscter size. The thruster fires continuously from its 
initiation until the end of the test period. The thruster is 
ieeely stoOo+small. A five newton thruster, 1.e., Figure 2, 
seems to be a viable engine size, although, without more 
PeeemmatliOneelt is difficult to, determine positively. 

Relatively large thrusters burn for shorter periods 
of time than smaller thrusters to achieve a common result, but 
each burn expends more fuel. Smaller thrusters, on the other 
hand, expend less fuel for a given burn time, yet must burn 
longer to achieve the same results. The fundamental result of 


Mmpemrecsct 1S that thruster size is a tradeoff variable that is 
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to be used in conjunction with other factors to achieve 
desired results, such as minimization of thruster burner 
during orbital maintenance. 

3. Final Tests 

Previous sections have shown that within a narrow 
radial bandwidth, a model using thrust equal to dGpagiiae 
superior to one using fixed thrust-angle reboost techniques. 
The question arises as to whether these results will remain 
valid for larger bandwidths. 

As described earlier, the control law commands the 
thrusters to fire when the spacecraft orbit has decayed a 
certain distance below the reference orbit. Results are 
examined for cases where the spacecraft is allowed to descend 
20 km, 30 km and 40 km below R,. Thrust is fixed at S0Gam 
while specific impulses vary between 300 sec, 500 sec and 2000 
sec. As for the case of a 3 km bandwidth, a thrust angle of 
70° maintained the spacecraft within the desired bandwidth and 
was therefore used for all cases described below. The tests 
are performed over a period of 100 orbits. In order to reduce 
computer time, atmospheric density was once again increased to 
a value of 1 x 10°? kg/m’. 

a. Case 1: Thrusters Fire 20 km Below R, 

Rather than choosing a specific bandwidth and 
adjusting the control law to achieve it, the spacecraft’s 


orbit was allowed to decay a certain distance prior to 
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activation of the control thrusters and the resulting 
bandwidth measured. This provided expediency since the control 
law determined the bandwidth rather than having to be adjusted 
to achieve it. The results are the same in either case so no 
@eeuracy 1S lost with this method. The figures in Appendix I 
illustrate the results of this case. 

MITE WwwenMer Elle = SpaCeecralt S Orbit” toe decay 20° km 
PremeretOo activation of the control jets”provided a radial band 
of 55 km. Plots of expended propellant mass versus orbit are 


provided for the three test cases ial 300 sec, 500 sec and 


sp 
2000 sec). AS in previous cases, thrust equal to drag provides 
a straight line while reboost is actually a series of steps. 
Fach step is a cycle wherein the spacecraft descends 20 km at 
which time reboost occurs (vertical portion of plot), after 
deactivation of the thrusters, the spacecraft once again 
descends to the point where reboost reoccurs. This is 
mmemmeared by the horizontal portion of the plot Since no fuel 
is being expended during this portion of the trajectory. As 
before, the reboost maneuver puts the spacecraft into a 
Seeomely eccentric orbit (on the order of 10°). This accounts 
Meme Periodic motion and high number of thruster firings 
indicated by the mass plots in Figures 2 through 4. As 


BePectea, increasing specific impulse reduces the amount of 


propellant expended over the test cycle. 
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b. Case 2: Thrusters Fire 30 km Below R, 

Requiring the control law to fire 30 km below the 
reference orbit provides a radial band of 78 km. Figure 5 is 
an illustration of spacecraft radius over time. Because AR is 
so large, the reboost maneuver must occur twice before the 
Spacecraft settles into a periodic trajectory that carries it 
the full width of the radial band. Eccentricity (oie 
osculating orbit, however, remains on the order of 10) 930@amee 
in its periodic trajectory, results are very similar to those 
of case 1. Figures 6 through 8 illustrate propellant mass 
consumption over the test period. The first two thruster 
firings are clearly evident. Once in its periodic trajectory, 
the mass plots are very similar to those of case 1 and occur 
for the same reasons. 

Comparisons of the mass plots of case 2 to case l 
reveal interesting results. Percent difference comparisons of 
case 2 to case l, for identical specific impulses, provides an 
indication of established trends from which inferences of 
future results might be made. 

We see that for case 1, for each variation of 
specific impulse, reboosting the Spacecraft requires 482 
percent more fuel than the use of thrust equal to drag 
techniques. Case 2 reboost maneuvers require 578 percent more 


fuel than thrust-cancel-drag maneuvers. 
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c. Case 3: Thrusters Fire 40 km Below R, 

Allowing the spacecraft’s orbit to decay 40 km 
Prior to initiating reboost sequences provides a radial band 
of 100 km in which the satellite trajectory is maintained. 
Figure 9 illustrates radial deviation of the spacecraft. The 
increase in AR coupled with a fixed thrust requires the 
Satellite to perform three reboost maneuvers before the 
GPeamaliar periodic trajectory is achieved. Thruster firing is 
clearly evident in the first two incidences, as is the ensuing 
eraioat a] decay of the resulting (ce lege itedligy,) eccentric 
mrajectories. 

Calculating the mass expenditure percentages 
reveals that reboosting the satellite requires 481 percent 
more fuel than does setting thrust equal to drag. Although 
thrust equal to drag is still Superior, the difference between 
the two is decreasing. To further test this result, 
percentages were calculated for points 43 orbits and 97 orbits 
into the test period. All values were less than corresponding 
values calculated in case 2. While these results do not 
provide concluSive evidence, we may conjecture that a trend is 
developing, indicating that at some point results from 
periodic orbital maintenance will equal those from forced 
Keplarian motion, or as the theory predicts, the plots will 


reverse and periodic thrusting will become superior. 
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C. SUMMARY OF RESULTS 

Initial testing, utilizing a constant atmospheric demeua 
model provided a baseline against which, further testingue@ewess 
be compared. A thrust angle of 70° was found to produce the 
desired results. Fixing the thrust vector at this angle 
maintained spacecraft orbital deviation within a three 
kilometer band nearly indefinitely. Upon determining this 
angle, the computer program discarded the constant density 
model and incorporated an exponential atmospheric density 
model. 

This Simulation was then compared to a model in which 
thrust canceled drag. According to the optimal control theory 
developed in Chapter III, the simulation (reboost model) 
should prove superior to a thrust-cancel-drag model (relative 
to the problem of fuel-minimization). In fact, “rege 
required significantly more propellant to maintain the 
satellite orbit within the three kilometer band than did 
thrust equal to drag. To test the robustness of these result, 
specific impulse was varied between 200 and 2000 sec and 
thrust was varied between 1 and’ 35 N. The results remained 
unchanged. To further test the results, the simulation was 
tested over much wider radial bands. The results still proved 
thrust-cancel-drag trajectories superior to reboost@iiga =a 


although the difference in efficiency seemed to decrease. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


As stated previously, optimal control theory states that 
orbital maintenance using a technique in which thrust cancels 
drag, is not optimal. This means that some scheme using 
periodic thrusting must then be optimal. Through the 
complicated techniques of optimal control theory, a thrust 
vectoring scheme is shown to indeed be the optimum. If the 
thrust vector always points along the primer vector, the 
trajectory is optimal. A sub-optimal scheme using fixed-angle 
thrusting and parametric variation is presented here as a 
Simplified method of determining the optimality of orbital 
maintenance. 

In each series of tests, minimization of propellent mass 
using fixed-angle thrusting has proven to be inferior to that 
in which thrust is set equal to drag. At first glance these 
results appear to contradict the theory developed by Ross and 
Melton [{Ref. 3]. For small perturbations forced-Keplarian 
motion proved to be superior to periodic fixed-angle 
thrusting. As the perturbations increased (indicated by the 
increasing size of the radial band), it would seem reasonable 
to expect that the difference in fuel consumption between 
these two techniques would increase. However, the results of 
tests described in section A.4 of Chapter V reveal that for 


large radial bands, the percent difference in propellant mass 


33 


expended between methods of orbital maintenance using periodic 
thrusting and forced Keplarian motion, actually appeareuaae 
decrease. Based on these results, we may conjecture that the 
percent difference between the two methods tested here will 
continue to decrease until periodic thrusting yields results 
Superior to those for forced Keplarian motion. Further testing 
is required before the analytical theory proposed in [Ref. 3] 
may be conclusively verified. 

The problem as presented here is that of fuel-minimization 
during orbital maintenance. Solving the Mayer optimality 
problem derived in Ross [Ref. 3] yields the primer vector. 
This is a very cumbersome method requiring solution of a two- 
point boundary value problem. If periodic thrusting is done 
along the primer vector, fuel will be optimized. This thesis 
has proposed a simpler method using the energy balance of the 
satellite to achieve similar results. Results, however, 
indicate that for small perturbations, forced-Keplarian motion 
will provide the best results. 

Ultimately, the goal of this thesis is to provide a method 
of fuel-minimization that is practical and may be applied to 
existing systems. Propulsion systems utilizing vectored thrust 
are highly complex and have a correspondingly higher chance of 
failure. Additionally, the extreme size of the perturbations 
required before periodic orbital maintenance would become more 


economical than forced-Keplarian motion is impraeereaie 
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Based on these conclusions, and with the added knowledge 
that a variable-thrust propulsion system capable of operating 
continuously over the lifetime of a satellite may also be 
impractical, it is recommended that further testing of 
periodic fixed-angle transverse thrusting schemes for small 
perturbations be accomplished. It is recommended that the 
Mayer optimality problem described in Ross [Ref. 3] be solved 
and the primer vector found. The results should then be 
compared to those described in this thesis to determine the 
actual amount of savings acquired through optimization. It is 
possible that the amount of fuel saved may not warrant the 
Seenmeor building a propulsion system capable of periodic 
Peete thrusting. Additionally, a comparison of orbital 
maintenance techniques presently in use, with results found in 
this thesis, should be accomplished to determine their 


relative efficiency. 


SD 


ORe Oa] OOO O00 00 0-000 00.0 OOM O@ O80 O2O 0807010 Gl ten ene ee 


APPENDIX A 


OBJECTIVE: DETERMINATION OF FIXED THRUST ANGLE TO MAINTAG: 
ORBITAL DEVIATION WITHIN A PREDETERMINED 
BANDWIDTH. 


VARIABLE DEFINITIONS: 


X (1) 
X (2) 
X (3) 
x (4) 
XDOT (1) 
XDOT (2) 


RADIUS (METERS) 

RADIAL VELOCITY (METERS PER SECOND) 
THETA (RADIANS) 

ANGULAR VELOCITY (RADS PER SECOND) 
TIME DERIVATIVE OF X(1) 
TIME DERIVATIVE OF X(2) 
‘DOT (2) TIME DERIVATIVE OF X(3) 
XDOT (4) TIME DERIVATIVE OF X(4) 

RO = REFERENCE ORBIT 

D = DRAG (N) 

EO = SPECIFIC ENERGY OF REFERENCE ORBIT 
E = SPECIFIC ENERGY 

MU = GRAVITATIONAL CONSTANT 

M = SPACECRAFT MASS (KG) 

GAMMAR = FLIGHT PATH ANGLE (RADIANS) 
GAMMAD = FLIGHT PATH ANGLE (DEGREES) 
TH = THRUST tN) 

TMAX = BLOWDOWN (MAXIMUM) THRUST 
ALPHAR = THRUST ANGLE (RADIANS) 

CD = COEFFICIENT OF DRAG 

RHOO = REFERENCE ATMOSPHERIC DENSITY 
RHO = CALCULATED ATMOSPHERIC DENSITY 
SPI = SPECIFIC IMPULSE 

V = VELOCITY 

SREF = REFERENCE SURFACE AREA 

GO = GRAVITATIONAL ACCELERATION 

H = INCREMENT OF TIME (STEP SIZE) 

PTI = PRINT TIME INTERVAL (STEP SIZE) 
T = BEGIN TIME 

FINAL (END) TIME 

a = SEMI-MAJOR AXIS 


Hh WM 
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OPC ©) 


om 
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n= ECCENTRICITY 


START PROGRAM 
PROGRAM ACTORB 


VARIABLE DECLARATION 
iMPbrcia REAL«=8 (A—-H, Lb-Z) 
DIMENSION X(4),XDOT (4) 


DEFINITION OF CONSTANTS 
Pos 1415926535 9D+0 
GO=9 .806D+0 
T1=0D+0 
N=4 


fp ee CC CC CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCE 


S 
ec 


Q 


MAIN PROGRAM 


Seen (10,FILE=’ INIT’ ,STATUS=' OLD’ ) 
Seen (11, FILE=’ OUT’ , STATUS=’ NEW’ ) 
OPEN (12,FILE=’ ORBEL’ , STATUS=" NEW’ ) 
OPEN (13, FILE=’ RAT’ , STATUS=' NEW’ ) 


READ (10,1)R),V0,M, TMAX, T, TF,H, PTI, CD, MU, RHOO, SREF, SPI 
memmAT (13(/,21X,D13.7)) 


PRINT*,’ENTER ALPHA’ 
READ*, DEG 
ALPHAR=DEG*PI/180.0D+0 


INDEX=0 

KOUNT=1 

X (1) =RO 

X (2) =O0D+0 

X (3) =0D+0 

X (4) =1.1673449D-3 
EQ=M* ( (VO*V0) /2-MU/RO) 


ioe, (| yeep’ TIMES ORBITS RADIUS VELOCITY ALPHA 


* ANGLM ENERGY TMAX’ 

eee l1, *),‘ (sec) (km) (km/sec) (deg) 

xi 

WRITE (12,*),° i; a e APOGEE PERIGEE 


ai 


O.@ eS 0) 2 Oxaere@ 


OG. © 


OQ GO 30-6) OO) 


NO 


QD & Ww 


*PERLOD 


WRITE (13707 TIME ORBITS DRAG TH 
VSS S) GAMMA’ 
CALCULATIONS 


CALL DRAG (SREF, CD, X, RU, RHOU, D7 aa 

CALL ONOFF (RO, X, TMAX, SPI,G0, TH, EO, GAMMAR, ALPHAR, V, D) 
CALL DIFFEQ(X, XDOT, MU, DM, MAX ee 

CALL RK4(T; xX, XDOT, NZ, INDE.) 


IF (INDEX .NE. 0) GO 7@ i390 


UNIT CONVERSIONS 


R=X(1)/1000 

SPEED=( ( (X (2) *X (2) ) + (2 (1) *X (4) 42 0 oe a 
V=SPEED*1000 

M=M- ((TH*H) / (SPI*G0O) ) 
MASS=20000-M 
GAMMAR=ATAN (X (2) / (X (1) *K (4) )) 
GAMMAD=GAMMAR*180.0D+0/PI 
ANGM= (X (1) *V*COS (GAMMAR) ) 
ALPHA=ALPHAR*180.0D+0/PI 
ENERGY= ( (V*V) /2) — (MU/X (1) ) 
ORBITS=T/5362. 458 


CALL ORBDAT (ENERGY, ANGM, PI,a,e, APOGEE, PERIGE, PERIOD) 


OUTPUT 

IF (KOUNT .LT. DNINT(PTI/H)) GO) [Cmz00 

WRITE (11,2)T, ORBITS, R, SPEED, ALPHA, ANGM, ENERGY, 1a 
WRITE (12,3) T,a,e@, APOGEE, PERIGE, PERIOD 


WRITE (13,4)T,ORBITS,D, TH, MASS, GAMMAD 


FORMAT (2X, F7.0,1X,F5.2,3%,F8.372%, Pome aor 


*F14.2,2%,F12. 2,25 oso 


FORMAT (2X, F7.0,1X,F10.3,2X,F4.3, 1%,F10.3, 1 3eb 102s, 
FORMAT (2X, F7.0, 2X, F5.2, 2X, F12.9, 1X, F5.0,0%) F153) 
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200 


& 


Cc 
c 


KOUNT=0 
KOUNT=KOUNT?+1 
twee. GE. T) GO TO 100 


END 


Gummer CCC CCCCCCCCCCececccecceeececececceeccececeececcccccceeéec 


c 


C 
c 


SUBROUTINE DRAG(SREF,CD,xX,R0, RHOO,D,T,V) 
iP LIciT REAL*8 (A-H, L-Z) 
DIMENSION X(4),XDOT (4) 


Boe (2) *X (2) + (X(1)*X(4)) **2)**0.5 
BETA=2.12D-05 
RHO=RHOO*EXP (-BETA* (X (1) -RO) ) 

0. 5D0*RHO*CD* SREF*V*V 


RETURN 
END 


pee eee eo COCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 


g 


PWUBROUTINE DIFFEQ (X, XDOT,MU,D,M, TMAX, ALPHAR, T, TH) 
Perici’ REAL*& (A-H, L-Z) 
DIMENSION X(4),XDOT (4) 


A=X (1) *X (4) *X (4) 

meu / (X(1) *X(1)) 

Meee 1) = (XZ) / (1 (X(2) *X (2) ) +(X(1) *X (4) ) **2) **0.5)) 
E= (TH/M) *SIN (ALPHAR) 


P=2*X (4) *X (2) /X (1) 
eee cs (1) *M) ) * ( (X(1) *X(4)) 7 (CCX (2) *X (2) ) + (% (1) *X (4) ) **2) **0.5)) 
S=(TH/ (X (1) *M) ) *COS (ALPHAR) 


XDOT (1) =X (2) 
meet 2) ~A-B-C+E 
XDOT (3) =X (4) 
mDOT (4) =-P-Ots 


Bo 


o9 


C 
e 


RETURN 
END 


COCCCCCCCCCCCCCCCCCCCCCCCCCC COCO COC CC CC e ee e eee CCC eee ee ee 
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SUBROUTINE RK4 (T,X, XDOT,N,H, INDEX) 
IMPLICIT REAL*8 (A-H, L-Z) 

INTEGER INDEX, I 

DIMENSION X(4),XDOT (4), SAVED (4) , SAVEX (4) 


INDEX=INDEX+1 

GO-TO, 1,272,449) I hbex 

DO 10 I=1,N 

SAVEX (I) =X (I) 
SAVED (I) =XDOT (I) 

X (I) =SAVEX (I) +. 5D0*H*XDOT (I) 
T=T+.5D0*H 

RETURN 


DOeZ 0 — nN 

SAVED (1) =SAVED (I) +2 .D0*XDOT (1) 
X (1) =SAVEX (1). oD07 eben 
RETURN 


DO 30 I=1,N 

SAVED (I) =SAVED (I) +2 .D0*XDOT (I) 
X (I) =SAVEX (I) +H *XDOT (I) 

T=T+ .5D0*H 

RETURN 


DO 40 I=1,N 

X (1) =SAVEX (1) +H/6.D0* (SAVED (UI) bon Gi) 
INDEX=0 

RETURN 

END 


cecceececceccccccccccecccccccc CCCCC CCC EEC CC CCC ee Cee ee eC ee ee eee 
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SUBROUTINE ORBDAT (ENERGY, ANGM, PI,a,e, APOGEE, PERIGE, PERIOD) 
IMPLICIT REAL*8 (A-H, L-Z) 


e= (1+ (2* ENERGY *ANGM*ANGM/ (3.98601208133D+14* 


40 


mo 90001208133D+14)))**0.5 
a=(-3.986D+14/ (2*ENERGY) ) 
APOGEE=a* (1+e) 
PERIGE=a* (1-e) 
PERIOD=((2*PI)/(3.986D+05 


/1000 


ASO: 5) east ae) 


C 
c 
RETURN 
END 
C 
C 
Memes Se CCCCCCCCCCCCCCECCCCccceccccccccccccecececececcccccecececcc 
C 
SUBROUTINE ONOFF (RO, xX, TMAX, SPI,G0O, TH, E0, GAMMAR, ALPHAR, D) 
IMPLICIT REAL*8 (A-H, L-Z) 
DIMENSION X(4) ,XDOT(4) 
é 
Cc 
DELTAR=R0-1000 
mec) ~ X42) )+(X(1) *X(4))**2)**0.5 
MU=3 ..986012D+14 
M=20000 
B=M* ( (V*V) /2— (MU/X(1))) 
Cc 
@ PRINT*, ‘MU=',MU 
Cc PRINT*, 'E0=',E0 
c 
IF (TH .EQ. TMAX) GO TO 99 
c 
IF (X(1) .LE. DELTAR) THEN 
C 
imine LT. EO) THEN 
TH=TMAX 
ELSE 
TH=0D+0 
ENDIF 
ELSE 
TH=0D+0 
ENDIF 
c 
c 
99 IF (E .GE. EQ) THEN 
TH=0D+0 
ENDIF 
c 
c 
100 RETURN 
END 


41 


— 


PECECCCOCCECECE 








CECCCECCECECeet 





APPENDIX B 


Z eanbtg 


u@WO 


Oe gt Qi 2b Qi Gi wl Cb Zt bt Qi 





on ent #8 NHN = © 


€ eanbtg 








(uoy) 
Smudve 


- 


T eanbtg 


LUGO 


, b @b 2b Qi Qk OL Gb Jb bb Qs 


(n) 
LSOMHL 





(ovr O14 x) 
AOU3N3 


43 


g ernbdty 








LIGHO 
“ 
80 
(uo) 
Snigi fd 
ae 
ase 
oe 
G eanbdty p ernbty 
gO 1640 
° be Gb @t Zt Qt QL OL EL Zt dt QL os 
: <Qe0nons- 
g 
¢ <A oTDIOE- 
K! <(ngw0DS- 
: (N) @ur or ») 
» tenuHL COLO ARENA 
‘ <ogsg00s- 
“ 
<ONegOK- 
<URpeNIS- 
OamBO- 


Q=-® 
€. uy (g-)a1 «= 4 
NOt @ aL 





44 


APPENDIX C 


€ eanbtg 


41HO 


oe Qi gt 4) Mi Gi wi Gi Zi tt Qi 


Z eanbrtg 


OF gi Gi 2) Qi Gi pl Gt Zi bb OI 





+0 
b 
& 
& 
? 
e 
cs] 
ra 
8 


ol 
bb 
a 





bigatecacaaae 
? 


T 9anbta 


LUGO 


‘ b @) 48 @) Sb wi Gb at bb Gb 


(N) 
Isnt 





our o1 9) 
AQt’3aNna 


45 


g eanbdtg 


LiGRsO 


oS @i Qt 4b Gt Gi wt Gt zi bi Gs 





(nw) 
LSNWHl 


eo ~ @a Ff #8 WwW = © 





p eznbta 


igyuO 


OS Gi Qs 4k Qt 9b wi Ct Ai 44 Qs 





vr o1 4 


46 


APPENDIX D 


LIGkiO 


' b qb 4b Qt Qt wi Gb ai bb QL 


€ ernbtg 


ugiO 


OZ qi Qi 4) Qi Gi wi Et Zi bs Qh 


¢ eanbtg 





(Nn) 
ISNiHL 


, 
1 6 





T ernbtg 


LIGuO 


ee @i qi 4k Qt Gb ot Gb ab bt Qh 


t 
Lo 5 


(ur o1 0 
AOW3NA3 





47 


G eaznbdtg 


LugidO 


DS qt Qi 2b Qi Qi pi Eb Ai bb Qh 





(vr o1 =) 


48 


g ernbTg 


LGdO 


oS @t @t 2b Qt Gi wi EL 2b bi Qi ¢ 


gO 


OS @t Ob 4b Qt Qs wt Gh at bb Ob 


Hi 


+ i ] 7) " ¢ ‘ i ! r 
‘ 4 " 


(n) 
oe: .JSNYHL 








6 ernbtg 


i 
7 18 


(un) 
ase SNiVU 


Gur 01 
ADHSN2 





PEEEEEEE 


i= ® 
6. uy (g-)a1 = d 
NOOS © mew) 


49 


APPENDIX E 


z eanbtg 


uO 


' L gb 2b @t Qb ps Cb 2b Lt Qh 


1, ua 


=» 


Gu¥O (1 4-)a0C 168 = 9 
NGS = 20Bu | 








im rT 8 # 


WN) 
LsnditHt 





' yi) 
(uo4) 
ose sfu0Vi 
— 
: 
ae ome 
6.uNlm (1 1-)30615 8 = ¢ 
NSS = EW 
{T eazanbdtgy 
JuGuO 


oS Gt Gi 24t Qi Gt OL Ci Gi Lb Qu 


mwa B 
§ uf (11-)30616 9 = 4 





N9Z = Meu, 


@vr 01 *) 
ADW3IN 


518 


OF qi Gt 2) Ot Gt pi Gi Zi 4b Qt 


g erndtg 


LUGO 


Rucee ee 


g orznbtg 


LugdO 


OS Gi Qt 4) Qi Bt OL Gt Zt tb GQ 


1H ek 


©. uy (1 4-)506 169 = 4 
NOZ = Mow ft 





* me 

agg 

aus 

9 2688 

(n) 

ess 2 SNK 

2G 

~~) 

eo) 

cma 

DL= 8 
yp eanbdtg 
LIguO 
OZ Gt Ot 4b Qt Qt wi Ct Zi bt Qh p 
<oScos- 
<8.a008- 
<0012008- 
Ose 00s 
«ops «= Bur 01 4 
ADU3Na 





DL=® 
g. wy (1 4-)a0E 16 B= ¢ 
NGZ @ MULL 





ot 





6 eanbdta 





LiguO 
' 6 @b 2b QL Qt OL GL at bt QL ? aces 
e008 
; ’ 
‘ ‘Lo ’ isg8 
@ 4888 
(uo) 
ee snivy 
Tee 
eS 
ij ; ry : 
t uy Oe 
Cae 
g eanitg 4 eaznbtg 
ugiiO LuQuO 
4 4 @b Zk @b Qt ph © 2k bt Qt By i Gt GL 24k Ot Ob PL GL Gb bt Qt 
z 
v 
8 
t 
ot 
et (N) 
LSNHL (Bur o1 
vt AD3N3 
91 
ot 
3" 


SZ a8 
¢. uy (1b-)S0c 16g = 9 
NGZ © XBW) 





a2 


APPENDIX F 


p ernbtg € eirnbdta 








Lugo LuguO 
. BG) Gt Li Qi Qi ws Gi ab 4b a pes 
¢ 
a ! ~~ 
° i ‘ ft r ‘ 
ab iil ’ } bite) 
vs ru 
es thr) (usd) 
2 INWi1n3d0tid 
on 
oc 
* — 
= ' 7 T : 
ec Thee epee 
‘He eE SS veces 
w 
0" ad — 
3 - = 2 =» 
we soe LE Bi OF guy (1 4-)30€ 16 9 = OF 
NSS @ Mow 
Z eazanbtg T eanbtg 
LuO iauO 
‘ ‘| @i 21 Qi Ql pl €b 2) tb Gl a Gg go OF Gi Gt 2b @t 9b wi Cb Ab bb OE ’ 
| - ee 
o <0a2008- 
mane 
m (vr 01 *) 
‘ 
e (nN) aap ADU3NG 
o)) | ASMWHL <09ga005- 
me <c0aa00s- 
e 
| a 
| on <openos- 
| ez <Ostalos 
Ps <O0ca00s 


2-8 


Cuddy (1 5-)a0c 16 8 = OF 
NGZ = XU | 






335 






AL = L 
§.uydy (1 b-)S0E 6 @ - OF 
NGZ = Ow) 


g eanbdta 


LigksO 


Oe Oi Bs 2b Ot Gt ph Ci Zi bb Gt 


Oie-® 


¢. ut (1 1-)30c 168 ~ 0° 





NOZ = ew) 


(Ox) 
SSVN 
LNVTISd0td 
to 

3 

v 

8 

8 

OL 

ze WN 

ASNBHL 

v 

Py 

e 

oz 

wa 

ww 

w 





L eanbtga 

LIGksO 
os Gt Gt 4b Qt Gb wt Gi at 44 QE ? ease 
4 toc 
‘ i) 
el 
ill 
— 
Smal 

y * : 
—a- 
ai 

G exrnbdtg 

16u0 
MS Qt Qt 4b Qi Sb wi Gb Zt bE QE + <osea0cs- 
bit tty 
<~oe.az008- 
<0.2008- 
ea 
<O0B 00s 
<osaa00s- 
«osm 
“eras 
<osecos- 
ey 6) yy 
ot 





Ai=*® 
€,uydy (11-)30€ 168 - Of 
NGZ = UL 








| 


ur 01 ») 
AOW3N3 


94 


TT eznbdtg 





4iguO 
‘ (} as at Pb Eb ab bb Ob a 
a 
ot vt 
‘ 4 
7 : "Ea ‘ 1308 
es 0.980 
(By) (un) 
= Scvn ie shucvu 
iW INV T13d0td 
ow eee 
oe | 
= io) 
'] 
oe 7 4 BY 7 nl 
) 
as (Cl 
Ilia BS 
& 4B (11-)306 169 = 07 
NGZ = KUL 
OT eaznbdtg 6 eanbdtg 
Liao 
OZ Qi @i 4} Ob GI ws GI Zh 4 OF a a. ' 
0 1S Gi Gi 4b Qh Qh wh GE aL 4s QS 
Zz —— 
y <opssu08- 
9 <0a.z008- 
” <0. a008 
ray <QuBI00E- ~<a 
ay 
rT coecu0s ABUAN3 
vi 
ds 
ge 
| | oz 
| - 
a 
w 


8d a 8 
€.uydy (14-)30C 168 = o4 
NG = Mew 


Gia B 






GUO (1 4-)306 168 = 07 
NSZ = MULL 


Dy) 


APPENDIX G 









386 (Op @ Cs 
¢ UME (1 4-)R0C 16 B= 


NGS = Ww) 


yp ernbtg 


gp 


© oS 
~ 


& 
pat 


29 


¥@s 8 Rk R 


¢.uiOs (41-)30616 8 = 04 
NGC = Mu) 





SBFREBSSSSSERRERSSCOS 





22 


o)5 


APPENDIX 8H 


















28 








oie) 


APPENDIX I 

















i 


“ii 


mu " 


i 


ie 


‘ 


ii! 


ii 


sae 


il 


° 
| 


us eazanbt 7 


ne ii 
ae 


He 


NE 7 


oe v¢ wp fo Of & © 











| 


h 








| 
| 


| 


i" 


Ie : 


j 


' 





li 


i 








Dye 





0001 


oop 


(Oy) 
SSvw 


4 


i 


me 


‘ 


i" 


i! 


i 


| 


i 


| 
| 


Hl 


Ry 


ini i 


i 


Ki 


vn 


ul 


' 


ih 


ne 





: 


| 


Wii 


; 





le 


| 





| 


i 








iy 


’ 





hi 





Se, 


ZT earnbta 





TT e2nbta 


6 eanbta 


LuguO 


oo oe ef eS Do eS e HY FE 


( 








A 


i / a 
ni % 
all 


634 


i i 
ii 


i ne 
OS 





EGER ER ERS 





(uo) 
snidvd 


60 


LIST OF REFERENCES 


Lawden, Derek F., Optimal Trajectories For Space 
Meavigation, Butterworths, London, 1963. 


Marec, J. P., Optimal Space Trajectories, Elsevier 
Scientific Publishing Company, Amsterdam, The Netherlands, 
1979. 


merece tf, M. anc R. G. Melton, Singular Arcs For Biunt 
Endoatmospheric Vehicles, Proceedings of the AIAA/AAS 
Astrodynamics Conference, Portland, Oregon, August 20-22, 
1990, AIAA Paper # 90-2974. 


Sm, Christina L., Orbits Conteining Arcs of Minimum 
Altitude Variations, Masters Thesis, Air Force Institute 
of Technology, Wright-Patterson Air Force Base, Ohio, June 
1989. 


Eisner, A., Drag Estimation and Satellite Orbit 
Determination, Applied Physics Laboratory, John Hopkins 
Mmrversity, March 1980. 


Kaplan, Marshall H., Modern Spacecraft Dynamics & Control, 
wiemm Wiley & Sons, 1976. 


U.S. Standard Atmosphere 1962, p. 57, United States 
Government Printing Office, 1962. 


61 


INITIAL DISTRIBUTION LIST 


Defense Technical Information Cemeer 
Cameron Station 
Alexandria, VA 22304-6145 


Library, €¢éde 2 
Naval Postgraduate School 
Menterey, CA "33943 [5002 


I.M. Ross, Code AA/Ro 
Naval Postgraduate School 
Monterey, CA 93943-5000 


David D. Cleary, Code Ph/Cl 
Naval postgraduate Scheel 
Monterey, CA 93943-5000 


David Pauls 


VAQ-129 
NAS Whidbey Island, WA 98278 


62 


No. of Copies 


Z 











DURE mise tee ene 


NAVA ‘HOOL 
MON?T®. 101 








1 te 
“Or gels ie in’ erebae 
ah fede ots 










re 





DUDLEY ITT nV “ a 


oe es Hl tl HH MUA IK I i Wieey. Sree oe a 


2768 00307708 2 epee ee ts 














- £ penton 
942 ignite bp Anthea dl & 
Seat a a 









a hy bar} 
Mi, bEmtive 489 0H 
Dhed date © phe ep 








one haeey eras Pall 6 4 
Ld ¢ 

ene oe 

obs tom ued? 

otade 


‘ ragete ate 
tat ted Soagises 

et a rh fie y we 

ee 


ni hignths 5 








































Ls fae: 
dat .aee! Fe i sat RP 
ae we Be 1 RNs Od ' ar ole nerds Loose fat o'e Sept 
Boe aaa ter AY et tate Teh) a Doe Y, ih poaes 
's steed > BPP ? eslitg tal w 
. BEES 
ae 
oa 


a ipapet ict 

rps sepsis Paears AD “Sr agra 
. . . 

a a thangs Oe oh as ei dyed teh . ha hg 

bf ‘ saline prs any * 

Witahalate?, ile Sree Dg 

aloud ’ 





See oe" at 
wk 























a Ty 

rors 
a Hei tty z 
5 ote bead Beate mr suc ede 
®, ie adnhe ty: 

4 Ped bs tnt 2 


Gee mde 
ved! eT en ies 
tated CA Mot Staal & Stee 

PA ee a dai 
towree Was forte tg! ‘. 























































































ee 1 8 


Pid Tels 
aytge fer, 
ea bY " 1 
eo Ney ek oe hon ines) ware . rage Ss 
1 . fogvane) ts et 1? ‘ ep F 
Pe . ‘ me ia 
' ‘ OO ee | we! ree a : 
Pate Ouagben ‘Pe ae a | ‘ = Ro 
. Peery 
Serie a al ae iSite ee | Pee : ma i. ae 
ern ¢ 8 t@ ’ ‘ ed@ntpes . . 
Ade pote OO a 8 ‘ ra oun . . ‘ 
‘ 
Oe acerca fa HL | . a 
ye Nae ae tek ald wee . oe 
' t efae | 
fy were Oy ’ hee x i 
‘ 
Sst aha oe . 
ee 
. . 
‘i Wt 1 . 
a ! be 
a teat pS . 
. € 
| at . , 
: : e ' ‘ 
ut pe Ph ri . » 8 
eo oh a a8 « ' ’ 
Lwmaettinnte & en 7 
dg Sant ge he 4 nat, i 1 . 
Hees ak Byes ray es . R . . 
Poe eet > 
wohseee sa tiead aap Pubng cocaine: 
Dh vary» 1 te Aas yaad wae bate 5 aa, «8 aveh ¥: 
fo et onde wok I 2 ° 1 1 
eoixee gfe i ie aay 
“as 
" 0 ay cv Lamar et eee nen oar SS Pact eer ; sae 
en ry One to a . 38 a, anes ' 
nA eon e e ‘ 


Bray ante Sf 
Meace be AAatet bs 
ode °ve se deer * ad fet plow 

- sb o* "a heel o> 
ve éal 















aaa 4.0 , “Wy nye ‘| 


pape Sour 
vel 2 fered, 


“a 
at behets Petri 
dV AA BRD eed! Fe heey 00 0%8 













pees | 
ari t ebad 








































ae awh? of: 03" ie baners 
he! ops) 
“a etass aoei Se i by 
Pune Veet 3 oe fips icrte im «tongs % 
Sh tte y are ry a2 he eT eS oe one 
RESET syenzerse nt wees pt ane 
I gh vice Ae 


Ty at 0 arei et of 









A amis deo be 16 bet 


ihe 













4 
R ees prrd act ' + al: . i 
Baten ive tb suahc 
wees é . 













Ars 





an ow *S* 
rah Bons ede 















gee 4 













ae 
c 2 
’ ‘, DSB et ‘Peete 4 
> 'e > whe ee 
ia Po% tet ae oh *. tiedusrt’3 te hiSeaabe teh Aye aly ‘ 
eres Lied be opel 23 Say hy 9 *s “dy tQue ae wre a 
243 6 ad, TARE. aPee Lives ,'*4 
wets peahyrew eget 
» Pah WLW 


‘ 
“stae ooo! yo <0 
' tes ab a's 
wed? . 


































© Sees 
aBe Saath: hed os 
7 
















Dory FIG 
Ie ue we a bee be? 
adh ate vedanta’ odd ares M mie Ve gut Sak oP 
F¥ su. eal. Red bphes dt oPy he 0 8: oy8. , Ss 
* pete > a ag . i ivy ism Sd ae'os 2 

i a; : sedan Prayte His teares fro i miton ry fy ¢ 
art eeei? ne oblige’ rAd Gea SO ep? : 
ler "t ogre : 


























a 

‘rest mec 
. SU . 
Piet g 










(i eee 
oe Ni rpeass 
aS 


te 









“9! 
yen Per ath 
s 





Ron » : Essh ts 


af oye 
ER i Peotasa + Asp Mhe sonra 
































¥! eRe ene Bre faen® ° 
s ei dette Me tay Oe isch 
Ay we en oe Soper reed ends ? 
re wv Beas Tepes 
athe ee A ony ea ~ Joye i? e 
on ee fe « * Noe hile e 
re 
° eos 
. o7 e« 
ot) 1 * 


Tae gee thee Ree 





ftltes 






ob gyigettene ie) 
al 5 43! 
Os ys hw 









iP re Rages 
1@ he vata, 
~f tobaled cher ® 
















5 fe Pace e 
s*, 459, Sb 
Saseeee 3° 







nee se Sn 
wrt * 

>, = Mine Seat bata 
mre Boar 1) pipes Lod. Yes. 44 
Smiryesth = oP! Aesctp preter 5 
staat ase pedi wont Be 


ao 




























Prete 
eds es? safe Aah 


t 


oe OM ? 
rOyss es, 
» 


Ty 


eet re wits 












t Mend eto 

















on 
at eli 











































pri Sentae, 
\é@ . 
H r Orin Fi 
a eee 5 38 nied BEANS cel spree g 
Poe! Lea '. ' 
Eee Stee Bah eee, Sale fia MSE teae OD 
: ? . 5 : . 
4 eS i .. adhe aed wy th % A. oe et we Sirs FS » 
we A rere, Trak te 2. a Etheh whey ' 





Toe LEAT 































































































































































































































































































































































































































































































ha Reto 7 apertures tare 
+ 8" ae ee o.. 8 ae Ob ww Sra cipae: Lay 
‘, fre EE ae Lal Lad i 3 
Cae ce} PE Aahs Fb Reh Bhs” ty. F ryt. 
3ee. Aue Fe RS rt. AS" : 
oH try Al a : 
a) : an > cere Hare os 
2 ated a! ° ae Te Cee re ; 
« Sere bive tT. pias, ¢ bos tee hed 
is tas ¥ BUEN ys BF ode nce. 
t Mantel iPttd % 
. Tyeidee> & OTF tt vt pee eS oD rested spree 
rs Wea a? eyls ae a PG e 
‘ atte ts Cre ost tte 
Paes ae Petter gd ee 4 
or) ‘ nee . 
‘ * » et x ba) . ' 
. « ae + 
Re ate e on aoe oa 8 
aha are Whee ew SAD Saget ee " . , 
errata feels 2 § a e. Bas 
f VLD oF Lap a Na eer Sev eeetare ss oa 8s - . soe oe tee 8 . . . 
b% eh ef my etsrery ye 4 : hy meg, © 8 Ope use 8 8 oo Os i one po eee as 
bya, may Ne he aey site eee £14. eh ies ee Es * Sak a th obs "s Beisesieboe = cies 4) Saege see owe ’ whe 
vay 5 "senmel! 33, tyi-on Sy Sarg “¥ nse tas mt ‘ : A : S pedneee 8 43 ers et eC Pe cn PT a ee . witin . 0! a 
net sare 13 4oM = vs Cas Sa) 34 ere oo Pea ty Be evar es er" Lh " ; - bd  @) . . Be a ae Y Pow 8 
mioed ys Ne cay yess ert Ava lg Stet ai , 834 2 a! ¢ ’ ‘ ; wd ’ * a eg 1 e a 
at es thik io cx tatate tobe ables ret * arect. oon’ s i¢ , ai P ze) : ; aa ; 
. safety Bas fy oy bey AESOP | athY * 33 MEeytm, oils t : 5 ' 
fcr ° Rad Las tt sib aye . ’ . e 
wag SE ws paage iC 77, NneP O Paspis, ee % 8 e 
Airey NUS ; ih st spans; Be j ais Mette - 
au ts es Soa giar ee? oyty® sy btearede ans wg 
fu dy i sare Bh reeh, ma 2 . - . 
siete “eS Mot . . a . 
Siratcns A: 1 : ’ 
Rte 2 At 
oonety® airrerey te oot hes . 
tat ve Amipht SEF % yd y oe So) 2 : 
Soret PERSPeil a* ey 22 Os ol i Le on 
any ofa: ae Us Ine oe ee 
a! Zolat re th : . - . 
a oe ed 1.” SF es-2e! ne . ae | 
trashy "6 ores A ptehar te ts a ee | - os 8 8 
She A GF ene > sd wile ri : Sa © 
& : ok F 
iE CGR RE: Se ut : 
Fst bel Wr dat Sana oo. @ ° . 
ts Aa ve ee ae a sper on ate - t owe 
om ge A oy Pa tee ge art afer nT fA poe Rote etra- . : 
ds be . yy ~a oa am ' ae 2? es 8 8 . + ° « 6 8 ° 
whe + oye “Fe cae . “ Bogen ue % . ts . e 
a0 5 “ay, 8 yee ot: bee se ! *- "8 . ot x . LE - . A x . 
af : ¢ ayte fe ea? Ss. wld & 2 ee oe * . a . eral eee . ° 
Brat Sed ERR poy oeges we tis ° 
ae fe i ere et Secekte F. Aled at f . 
ones pea ? Ry CS be +H ee 
ye, 2’d, . 
Cs Bean! {7 Lae ne . . 
; ath . 
rank nated a! a8 . 
AU ever es ne 2% ope . Sie 
Spr ae J yeryeuta ore 2e AC 2 wt, ew oe uray! = < © | 
r "erm ad » et ae te 2A oe : . hehe: ' 
af eveiz Rea nE ivearks aed a) . '. e« 8 
Ey sarees a raat 438% ne tege Cad. My otg ate see rh i: Yee’ ae : . : ! 
Fr MsdaFoe pres Rete ade * 489 84" ye! : ae ee ee ager eeete Al ' 
seal ests Laas gbaratel Pye oa rst 25 3 nyhe®, QTE 4, : om 
oe ee E 2 
mats ayy) +4 weil, tS aM 2 b5e%4 | eo 8 on 
aii se hes aa \ > ‘és re . 
ee any ger ee ae 
i t ar et ’ 
7 ss) ry Le : as i bie . ry eer ¢ , 
Pie 
F, 
>, a 5 
wR Sy es Nay betley Pes bee PPB OIR =» . 
Asx ye ey yyy en faiyre wy - 
prance Pe RERe BEN fede ny 
i paw pwede? . ae 
ey PANE os He it yes neh ate, vert : : 
Ly BA ahr ( st Fea of : outs . . ‘ 
‘ Valae y > ao . : 
| ae . 
wy be os “ st ». 
° oY . a Fae t 2 . nae . = 9 
t > : . . 

ac ee eg oe Ep accpe : : 

2 :* [eke ° ° . . . 
hh) eu . ea rs 
tt wap bay ° ats i . “ te 

? ae . ok ne ’ . s 
an , ate . 8 Tus a ” . » & : 
*f x \ ‘ ; P . ey . . 
oh ec *e . nf Bian pe? a . r) Be - a e 
rn ol, vom ¢ 3" ae “wy zo" a = Ps ae ote He fqts A zl a < ' 
S esd gis? » es of Ws 18kg s ie *, iyi, mer) Me a wig cet 89 ar ghete eS ‘ aes Rigs ‘ . 
} > vit, ree ee ry . i ht Hy Ags fe ; cae ie hide: Rie dete 4! * wt as »! e ! ae s ois = . 
he ane t he ah on ¥, ae ery a ain a 3 © ele ln ’ 
Ne sdess Liane Wy Fla (es 8 4 ergy tae be RE ots ah . ae > * Y aes A . 
ee iat tt Nout de Ae LM aise a 48 eS Ay ay fae * ‘ E: 
yea ry Ketan pt af An a es ate "pe sastes . ae - : 
‘ ¥ yA pity BEY br ma wee de 1 ng. ce ‘ wee o Pa ee 
> Sf fbi ity: he PL oe er “3 Velingounl. Cee eo © 8 ee 8 le . 5 ; 
vee tht gt ty dr, ae tins Te tokege pte an, ii “is ojkeer etete y Sei eeceea j . 
; ieee rate rae Roy yh fn °, 4 bees * 36. mates th. ota! oon ts tee 8 S ® . 8 ’ ’ 4 
++ i are Pieri dd bs qses ty? ° a §ye 2 » % t ‘ 
we a i Ne toliois al “SNS aoe $ LAS te. ed 7? tha “F -;! . a 2 . . 
er “EL. 180 Ff ‘ 
eh So 4 Ayr et Oa Seat nace : . ’ ’ 
ve? ® tw Ls 
. . . . eo ° 





